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Chapter 1 

Cinemática 



1.1 Introducción 

La cinemática es la parte de la física que se encarga de la descripción del 
movimiento. Esta descripción se realiza sin indagar en las causas que lo provocan 
(eso sería la dinámica). 

La primera pregunta es: ¿Qué necesitamos para poder describir el movimiento? 

• Un sistema de referencia , desde el que el observador toma sus medidas. 
Dos observadores en distintos sistemas de referencias dan descripciones 
del movimiento diferentes. 

Utilizamos dos conceptos fundamentales en cinemática: El tiempo y el espacio. 
Si tenemos, por ejemplo, un objeto viajando en el espacio, la descripción del 
movimiento es poder decir el cuándo y el dónde, o lo que es lo mismo, su posición 
como una función del tiempo. Para ello utilizamos patrones que nos permitan 
cuantificar estas magnitudes en unidades del Sistema Internacional (SI). 

• Tiempo 1 : La unidad en SI es el segundo. El patrón se obtiene a partir de 
una determinada frecuencia atómica: lsf> vq s = 9192631779 Hz 

• Espacio (o distancia) : La unidad en SI es el metro. El patrón se obtiene a 
partir de la velocidad de la luz en el vacío (c, que es constante) y el patrón 
de segundo: 1 m <h> c = 299792458 m/s 

Entonces, con una cinta métrica y un reloj sabemos el dónde y el cuándo — > A 
partir de ahí podemos describir el movimiento desde un Sistema de Referencia 
dado. 



1.2 Movimiento de una partícula 

La posición de una partícula, que en general es una función del tiempo, la 
denotamos con el vector r(t). En un sistema de ejes cartesianos, podemos 
descomponer el vector posición en 3 dimensiones en las direcciones X, Y, Z. 



1 E1 tiempo es absoluto. En física clásica el tiempo transcurre igual para todos los obser- 
vadores independientemente de su movimiento relativo 

7 
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Í(x(t), y(t), z{t)) — x(t)i + y(t)j + z(t)k si el movimiento es 3D 
(x(t),y(t)) = x(t)i + y(t)j si el movimiento es 2D 

x(t) si el movimiento es ID 

Observar que mientras r(í) es un vector, sus componentes x(t), y(t), z(t) 
son magnitudes escalares. Se relacionan con el vector r(í) mediante los vectores 
unitarios i = (1, 0, 0), j = (0, 1, 0) y k = (0, 0, 1). Estos vectores están definidos 
de modo que el producto escalar entre ellos es: 

i-j = i-k=j-k = 

i ■ i = j ■ j = k ■ k = 1 

además |i| = |j| = |/j| = 1, y son invariantes. Conocida la posición de la 
partícula, podemos deducir la velocidad (y) y aceleración (o). 



1.3 Concepto de velocidad 



¿Cómo cambia la posición r con el tiempo? Si después de un intervalo de 
tiempo Ai la partícula pasa de estar en la posición r (t) a situarse en r(t + Ai) = 
r(t) + Ar, es decir, ha experimentado un avance Ar, entonces su velocidad media 
ha sido de (v) = Af/At (La notación (•) quiere decir promedio de lo de dentro). 
Esta cantidad nos da información de cómo cambia la posición de la partícula 
en un intervalo de tiempo Ai. Pero aquí lo que nos interesa no es la velocidad 
media sino la instantánea, es decir, cómo cambia el vector r en cada instante de 
tiempo. Para ello tenemos que utilizar incrementos de tiempo infinitesimales; 
nos encontramos con lo que en matemáticas llamamos derivada: 



v(t) - 


,. f(í + Aí)- 

= lim t 

At^O Ai 


-r(í) 


dr{t) 
dt 


Ar 


\ f(t + Ai) 







w 
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El vector v(t) nos proporciona información sobre el movimiento inmediato: 

• Dirección y sentido del movimiento: Hacia donde se desplaza la partícula. 
Lo indica un vector unitario í„ = ||r, con |u„| = 1. 

• Módulo de la velocidad \v\: Cómo de rápido es el desplazamiento. 

Ojo! Tanto el módulo \v\ como la orientación u v cambian durante la trayectoria: 
En general dependen del tiempo t. Por lo tanto, v(t) nos dice cómo cambia la 
posición en cada instante de tiempo. Si conocemos f(t) podemos calcular la 
velocidad v(t). 

Por 'conocer f(t)' se quiere decir que se conocen (i) sus componentes carte- 
sianas o bien (ii) su módulo y orientación 

(i) Componentes Cartesianas: 

f(t) = x(t)í+ y(t)j + z(t)k =>■ 



= dr(t) = dx(t) u dy(t)j | dz(t) n = 
dt dt dt dt 

v{t) = v x (t)i + v y (t)j + v z (t)k 
donde se ha tenido en cuenta que los vectores unitarios son fijos, y por lo tanto: 



dt dt dt 

(ii) Módulo y orientación: 

Es decir, nuestra partícula está a una distancia del origen r(t) = |r(í)|, en la 
dirección y sentido que indica el vector unitario u r (t) = ^W. 

Entonces la velocidad la calculamos derivando con respecto al tiempo: 



di 


dj 


dk 




dt 


dt 


dt 




m 


= r 


(i) u r 


(*) 



v(t)= á fu r{ t) + r { t) M f 



Ahora el vector unitario NO es constante, depende del tiempo y por lo tanto su 
derivada debe aparecer en la expresión de v. 

1.4 Concepto de aceleración 

¿Cómo cambia la velociad con el tiempo? — > Es lo que la aceleración nos 
dice. Análogamente al caso de la velocidad (que es la derivada temporal de 
la posición), la aceleración la definimos como la derivada de la velocidad con 
respecto al tiempo: 



dv(t) d 2 r(t) 
a{t) = ^T = ^ 

La aceleración nos informa de cómo está cambiando la velocidad: tanto 
cambios en su módulo |í/(í)| (más rápido o más lento) como cambios en su 
dirección y sentido, es decir, en u v (t). 
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En componentes cartesianas: 



a(t) = 


d 2 x{t)^ d 2 y(t)^ d 2 z(t) f 
' dt 2 " ' dt 2 J ' dt 2 



= a x (t)i + a y (t)j + a z (t)k 



1.5 De la posición a la aceleración, y viceversa 

Ahora ya sabemos cómo, partiendo de una expresión general de la posición en 
función del tiempo r(t), podemos deducir tanto la velocidad como la aceleración 
en un instante de tiempo cualquiera, mediante derivación. 

r(t) A v(t) -£> a(t) 

El proceso inverso es el de integración en el tiempo: así podemos, a partir 
de la aceleración obtener la velocidad, y a partir de ésta, la posición. 

Si la función que conozco es la velocidad, entonces la posición se obtiene: 

-^ = v(t) => dr(t) = v(t)dt 

r(í) f t ,t 

df = v(t')dt' => f(t) = f + / v(t')dt' (1.1) 

f(t )=fb Jt a Jto 

Observar que, a la hora de integrar, aparecen constantes de integración 
que matemáticamente son libres. Sin embargo, para obtener una expresión 
que físicamente tenga sentido, hemos de utilizar las condiciones iniciales del 
movimiento: en la ecuación en la que obtenemos f(t) integrando la velocidad, 
hemos puesto como límites inferiores de las respectivas integrales el tiempo ¿o y 
la posición en dicho instante de tiempo r(t = t ) que denotamos como fo- En 
los límites superiores de integración ponemos la expresión a obtener r(t) y el 
tiempo t. 

Si el dato que conozco es la aceleración, entonces la velocidad se calcula así: 

-^ii = d(t) =*> dv(t) = a(t)dt 
dt 

v(t) ,t ft 

dv' = / a{t')dt' =¡>- v(t) = vq+ a(t')dt' (1.2) 

v(t a )=v Q Jt a Jt a 

Para establecer los límites de integración procedemos análogamente a como 
se hizo previamente. Para obtener la velocidad, integramos con límites inferiores 
el tiempo ¿o y la velocidad justo entonces vq, etc. . . 

Para deducir la posición a partir de la aceleración hemos de integrar dos ve- 
ces. O, dicho de otro modo, introducimos el resultado que acabamos de obtener 
para v(t) (ecuación (1.2)) en la ecuación (1.1), donde se saca la f(t): 



r{t) =r + i! v(t')dt' } _, . .. ^ , , '■' 

v{?) =vo + J to a(t")dt"j 



a(t")dt" 



dt' 
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1.6 Algunas relaciones útiles 

Hemos visto como a partir de una expresión de posición, velocidad, o aceleración 
en función del tiempo se pueden deducir las otras dos. Sin embargo en algunas 
ocasiones es posible que nos enfrentemos a problemas en los que, por ejemplo, la 
aceleración no nos es dada como función de t sino de la posición r. En tales casos 
no se puede hacer directamente una integral en el tiempo, porque no tenemos 
una expresión que dependa explícitamente del tiempo sino que lo que tenemos es 
a(r) ¿cómo se obtiene la velocidad o la posición? Recurrimos a las definiciones 
de aceleración y velocidad que ya hemos visto: 



1.6.1 Si conocemos a(f) 

4r = a =*• dv = adú _ ,_ _ 7 ,df _ ,„ 

at _ j- > v ■ dv = a ■ di— = a ■ dr =>■ 

v = % J dt 

pv(r) r-r 

v ■ dv ~ / d(r)dr. 

El termino v-dves un producto escalar entre dos vectores. Si las componentes 
de ambos son v = (v x , v y ,v z ) y dv — (dv Xl dv y , dv z ), entonces el producto escalar 
es: 

v ■ dv — v x dv x + Vydvy + v z dv z 

entonces su integral la calculamos como la suma de tres integrales separadas: 

v-dv = / {v x dv x + vydvy + v z dv z ) = - (vi + vi + vl)lj = - (\v(r)\ 2 - |u | 2 
volviendo a la expresión anterior: 

UW)\ 2 -Wo\ 2 )= [ r d(r)dr^> 

z Jra 



\v{r)\ = \ hl + 2 / d(r)df 



1.6.2 Si conocemos a(v) 

Mediante separación de variables podemos calcular v(t) a partir de la expresión 
de la aceleración como función de la velocidad. 



. , dv dv f* f v(t) dv 

a(v) = — =^> dt = ——— => / dt — 



dt «O) J t J Vo a(v) 



de esta manera podemos despejar v(t) y después mediante integración y derivación 
obtener, respectivamente r(t) y a(t). 



í l 


,r(í) 


dr 


dt = 


/ 




Jt 


L 


v(r) 
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1.6.3 Si conocemos v(r) 

El procedimiento es muy similar al del caso anterior: 

, dr dr 

v(r) = — ==>■ dt = —— -r = 
v ' dt v(r) 

así, obtenemos r(t) despejando. Derivando podemos conocer v(t), y, volviendo 
a derivar, la aceleración a(t). 

1.7 Casos sencillos 

Existen dos casos particulares de cinemática en los que todo se simplifica y 
podemos obtener algunas relaciones generales. Se trata del movimiento rec- 
tilíneo uniforme y del movimiento uniformemente acelerado. 

1.7.1 Movimiento Rectilíneo Uniforme (MRU) 

El movimiento rectilíneo uniforme se caracteriza por tener una velocidad cons- 
tante y, por tanto, aceleración cero. 

a = 0, v(t) = vq = cte 

Notar que todo movimiento con velocidad constante ha de ser rectilíneo por 
fuerza. Por tanto es conveniente escoger un sistema de referencia tal que uno de 
los ejes (cartesianos) coincida con la dirección del movimiento para simplificar 
las ecuaciones. 

Si escogemos el eje OX de modo que sea paralelo a vq, tenemos que, inte- 
grando la velocidad (ver ecuación (1.1)): 



x(t) = x +vo(t-to) 

1.7.2 Movimiento Uniformemente Acelerado (MU A) 

Acontece cuando la aceleración a que es sometido el sistema es constante: 
o(í) = ao = cíe. Un ejemplo clásico es el tiro parabólico, donde un proyectil es 
disparado con una velocidad inicial dada, formando un ángulo determinado con 
la horizontal. El proyectil se ve sometido a una aceleración constante, la de la 
gravedad (a = g = 9.8 m/s). 

En el caso general este movimiento no es rectilíneo. Sólo lo sería en el caso 
de velocidad inicial cero o paralela a la aceleración. 

La expresión para la velocidad se obtiene integrando la aceleración: 

dv(t) _ r {t) ,_ /■* 

= ao => / dv ~ aorfí =>■ 



dt j{T j to 



v(t) = v +a (t-t ) 



A su vez la posición se calcula integrando la ecuación de la velocidad: 

— — = v(t) = v + a (t-t ) =*> / df= / (v + a (t - t )) dt 

dt Jf Jto 



Fundamentos de Física I - Guillermo Sánchez 



f(t) — fo — / v Q dt + / á Q (t — t )dt = 

Jt„ Jto 

't 2 t 2 \ ft 2 f 2 

v (t-t )+a [ - - ^ J-a to(t-to) = v (t-to)+a i^-f-t a t + t 2 



f(t) = ñ o + v {t - í ) + -a (t - í ) 2 



Notar que habitualmcnte las componentes de la posición x(t), y(t), z(t) tienen 
distintos tipos de movimiento: En el tiro parabólico, la componente en la di- 
rección de la gravedad se mueve con MUA, mientras que la componente per- 
pendicular a la gravedad se mueve con MRU. 

1.8 Cinemática del Movimiento Vibratorio Armónico 

El movimiento armónico es aquel que viene descrito por la ecuación: 



-v .n df^t) 2 



con lo constante. 

1.8.1 Características del movimiento 

De la ecuación anterior se deduce que aceleración y vector posición: 

• son paralelos. 

• tienen sentidos opuestos. 

Como en general en cinemática, se extrae información física de las llamadas 
constantes iniciales: En el instante de tiempo t = ío, se conocen la posición 
r(ío) = ^o y I a velocidad v(to) = ub- Si son paralelos (fo || Vo)> entonces el 
movimiento es unidimiensional, y escogemos el eje OX. 

LL Ju n tí Jb O 

La solución de esta ecuación es una función seno o coseno: 

x(t) = Acos(ujt — 6) = ^4sin(cjí — a), I a = - 

Donde A > es una constante y 8 es la fase inicial o resultado de la ecuación en 
t = 0. Dada la naturaleza de las funciones seno y coseno, que están desfasadas 
por un ángulo 7r/2, se puede expresar x de las dos maneras, ajusfando el desfase 
adecuadamente. 

La velocidad y la aceleración las obtenemos derivando, como siempre: 

x(t) = Acos(u:t- 9) 

dx(t) 



di 



= v(t) = —Auj sin(o;í — 



Fundamentos de Física I - 



Guillermo Sánchez 



d 2 x{t) 
dt 2 



a(t) = -Alo 2 cos(wí - 0) = -lo 2 x(í) 



Posición, velocidad y aceleración son funciones periódicas de t y acotadas entre 
un valor máximo y mínimo, porque el coseno está acotado entre los valores 
— 1 < cos(u;í— 9) < +1 y porque existe un incremento de tiempo T que llamamos 
periodo tal que: 

x(t + T) =x{t) 

v(t + T) = v(t) 

a(t + T) = a(t) 

el periodo es T = — . 

A u se la denomina frecuencia angular. Conviene no confundirla con la 
frecuencia i> = A. 

Tanto la amplitud como la fase inicial se pueden obtener a partir de las 
condiciones iniciales (xo,Vq) en el tiempo t = 0. 



xo = Acos(—9) = Acos9 = 

vq = — Alo sin(— 6) = Alo siriO 
y como eos 2 9 + sin 2 9 = 1 , 



eos 9 = 



smf 



Xo 

A 



Alo 



(1.3) 
(1.4) 



A 2 A 2 c 



1 



A 



Por otro lado para obtener la fase 9 dividimos la ecuación (1.4), de vo por la 
ecuación (1.3), de xq- 



vq Alo sin 9 



= lo tan 9 



tan# = 



!•'() 



xo A eos 9 xqlo 

1.9 Movimiento en Coordenadas Polares 

Hasta ahora hemos visto sobre todo cómo describir el transporte desde un 
sustema de coordenadas cartesianas. 
Y 




* X 



En 2 dimensiones, las cordanadas cartesianas de P son: 

r = xi + yj 
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En polares, necesitamos un origen O y una recta 00'. La posición de P viene 
determinada por su distancia r — \r\ a O y en ángulo 9 de su vector posición 
con la recta 00', definido de manera que sea creciente en sentido antihorario. 

p 




O' 

Los vectores unitarios ahora son u r , en la dirección de r y ug, perpendicular 
al anterior y en el sentido de 9 creciente. 
La posición es: 

r = ru r 

El paso de cartesianas (x,y) a polares (r, 6), y viceversa, es sencillo: 

x = r eos 9 r = \Jx 2 + y 2 
y = r sin 9 9 = atan(j//x) 

Como ya se vio con anterioridad, en cartesianas el cálculo de v y a es sen- 
cillo: al derivar nos aprovechamos de la circunstancia de que los unitarios son 
constantes. 

f(t) = x(t)í + y(t)j 

m _ ^ + ^ 

¿Cómo se calculan velocidad y aceleración en polares? También hay que 
derivar, pero ahora hay que tener cuidado porque los vectores unitarios no son 
constantes, y por lo tanto sus derivadas no son nulas. 

f(í) = r(t)ü r (t), 

dv(t) d 2 r(t) _ drlt) düJt) drit) du r (t) d 2 ü r (t) 

d 2 r(t) _ , , n dr(t) du r (t) , d 2 ü r (t) ,„ „. 

La pregunta es ¿cuánto valen "^ ' y — ^} ? Podemos deducirlo a partir 
de las siguentcs figuras esquemáticas: 
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Üg (t + dt) 



u r (t + dt) 




En la figura anterior vemos que los vectores direccionales unitarios cambian 
cuando la partícula se mueve. A continuación se muestra un esquema con sólo 
los vectores unitarios y sus cambios: 

düg 




A partir de esta figura vemos cómo se deduce que 



\dü r 
dü r _L ü r 



\u r \d9 

$■ dü r || Üg 



dü r = d9üg 



dü r 



de . 



lü = Tt U0 



(1.7) 



Encontes ya sabemos cuánto vale la derivada del vector unitario Ü r . Además 
podemos calcular la segunda derivada: 



d 2 Ü r d 2 9 _ dd düg 
~d& = ~dp Ue + ~dt~dt 

¡penden de la derivada 
también podemos deducir cuánto vale esta última 

düg = —ddü r = 



(1.8) 



pero estas expresiones dependen de la derivada ^-. Según la figura anterior, 



de. 



(1.9) 



\dug = \ug\d9\ dug 

dug || -Ur J dt dt 

Entonces, el resultado de la expresión de düg/dt que hemos deducido en la 
expresión (1.9) lo podemos sustituir en la ecuación (1.8). Obtenemos: 

2 



d ¿ Ü r d ¿ 6 . 



di 2 



di 2 



ug 



(1.10) 



Ahora ya conocemos el valor de las derivadas primera y segunda del vector 
unitario ü r y las podemos utilizar para reemplazarlos en las ecuaciones (1.5) 
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y (1.6). Comenzamos sustituyendo la ecuación (1.7) en la expresión para la 
velocidad en polares (ecuación (1.5)): 



;(í) 



du 

, % ~5í ~dt 
_ dr(t) _ dr(i) 

~ dt ~~ dt 



d0 ÍÍ 
M U 



ü r (t)+r(t) d ^ 1 



til 



dr _ 



de 



dt Ur + r di Ue - 



A continuación reemplazamos en la ecuación para la aceleración en polares 
(ecuación (1.6)) las expresiones (1.7) y (1.10): 



dü r 
dt 



d-B-TT 

!i u e 



dt 2 



d ¿ 9r¡ (d8\ 2 „-t 



d{t) = 



- drí f> ü r (t) < Q dr (t) díl At) 



d 2 r _ dr dd _ 

a= dfi Ur + 2 dtdt ue + r 



dt dt 

d 2 í 



r(t) 



dt 



-ug 



d 2 u T {t) 

dt 2 



dB\\ 
dt Ur 



d 2 r(t) ídB 

dt 2 ~ r \~dt 



Recopilando, tenemos que la velocidad es: 



drdB 
> 

J dt dt 



dt 2 



a o 



_ dr „ 



de _ 

r dt us - 



y la aceleración es: 



d 2 r(t) (dd 

~dt 2 r [dt 



drdB 

2 

dt dt 



dt 2 



"o 



Recordar que la velocidad se define como la derivada de la posición (en 
metros) con respecto del tiempo. Análogamente, podemos definir la veloci- 
dad angular (que denotaremos como u¡) como derivada del ángulo respecto del 
tiempo. 

d6 

u = ~r 
dt 

y la aceleración angular (a), como segunda derivada de B con t. 

_ d 2 e _ du 
a ~lti 2 ~^t 



Entonces escribimos v y a así: 
dr 



di 



u r + rujug 



v r = dr /dt Componente radial 
vg = ruj Componente polar 



" =i dí- rw2 )^ + ( 2w S +rai ' 7í ' 



= 4vf — ríú 2 Componente radial 



ag — lio ^¡r + ra. Componente polar 



Como podemos ver son expresiones complicadas, más que en cartesianas. 
Sin embargo son útiles para describir órbitas, giros y otros tipos de movimiento. 
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1.10 Movimiento Circular: Descripción en po- 
lares 

Si tenemos una órbita circular, y situamos el origen de coordenadas en el centro 
de la circunferencia descrita, tenemos que la distancia r es constante. Esto hace 
que la descripción del movimiento utilizando coordenadas polares se simplifique 
mucho, porque ahora todas las derivadas de r con el tiempo desaparecen: En- 
tonces: 

_ _, í v r = Componente radial 

v = ruiue < „ , 

1 vg = v — rcu Componente polar 

Es decir, la velocidad sólo tiene componente polar, o tangencial a la trayectoria. 
No hay velocidad en la dirección radial, porque el radio es constante: la partícula 
cuyo movimiento estamos describiendo ni se aleja del centro ni se acerca a 
él. Conviene recordar la expresión que relaciona la velocidad lineal v con la 
velocidad angular lo: v = reo. 

La aceleración también toma una expresión mucho más sencilla que en el 
caso general: 

2 ^ I a r = — ru¡ 2 Componente radial 

a = -ruj u r + rau g < d \v\ „ ' , 

[cío = ra = -Tjr Componente polar 

Tenemos dos componentes de la aceleración: (i) la radial o aceleración centrípeta, 
dirigida hacia el centro: Esta aceleración no cambia el módulo de la velocidad 
sino su dirección. Según avanza la partícula, su vector velocidad va cambiando 
de dirección para que la trayectoria se ciña a la circunferencia de radio r. El 
módulo de la velocidad sí que puede cambiar, y esto es consecuencia de (ii) la 
aceleración tangencial, que tiene la orientación del vector u$. 

1.11 Casos sencillos de Movimiento Circular 

Los casos más sencillos de movimiento circular son análogos a los que hemos visto 
en la sección [1.7]. Tenemos el movimiento circular uniforme y el movimiento 
circular uniformemente acelerado. 

1.11.1 Movimiento circular uniforme 

No solo el radio (i?) es constante, también el módulo de la velocidad lineal. 

v = cte 

Entonces la velocidad angular también es constante: 

v 

lo = — = cte 

R 

y la aceleración tangencial o polar es cero: 

dv 

ae = Tt =0 
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Por lo tanto la aceleración total solo consta de componente centrípeta o radial. 



a = a r = —ru) u r 



Finalmente, para calcular la posición en polares, necesitamos el ángulo 9 
como función del tiempo. Lo obtenemos integrando la u>: 



de 



e(t) 



de 



0o 



Lüdt 



0(t)=0o+w(*-*o) 



El periodo o tiempo de revolución es el tiempo T que la partícula tarda en 
regresar a la misma posición, o, dicho de otro modo, en describir un arco de 
ángulo 2ir: 

e(t + T) = e{t) + 2tt =¡> 6(t + T) - e(t) = 2tt = ujT =^> 



1.11.2 Movimiento circular uniformemente acelerado 

Es decir, con aceleración angular a = a^ =cte. La frecuencia angular se obtiene 
por simple integración: 



a 



dio 
~di 



w(t) rt 

dui = aodt 

ljo Jta 



bj(t) = uj + a(t — í ) 



Integrando de nuevo deducimos la posición angular 6: 

6(t) r t 



U! 



de 

dt 



u(t)dt 



ta 



e(t) = e Q + wo(* - ío) + 2«(* - f o) 2 
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Chapter 2 

Dinámica 



2.1 Introducción 

La dinámica supone ir un paso más allá de lo visto hasta ahora, la cinemática. 
En este capítulo no nos quedaremos en la descripción del movimiento sino que 
nos ocuparemos de explicarlo, estableciendo una relación causa-efecto donde la 
causa son las fuerzas que actúan sobre un determinado sistema y el efecto es el 
movimiento. 



2.2 Vectores y productos 

Ya que vamos a trabajar con vectores recordamos: 

2.2.1 Producto por un escalar 

Si Ü es un vector y b un escalar, el producto lo escribimos v = bu. 

1. El resultado de multiplicar un vector por un escalar es un vector. 

2. El resultado es paralelo al vector multiplicado, u \\ v. 

3. Ambos tienen el mismo sentido si b > 0, y sentidos opuestos si b > 0. 

4. El módulo del nuevo vector v es el producto del módulo del vector que 
teníamos u por el valor absoluto del escalar b, es decir: \v\ = b\ü\. 

2.2.2 Producto escalar 

El producto escalar entre dos vectores A y B que forman un ángulo a entre ellos 
es: 

A-B = \Á\\B\cosa 

Por tanto, 

1. Si A y B son perpendiculares, A ■ B = 0. 

2. Si A y B son paralelos, Á ■ B = \Á\ ■ \B\. 

21 
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3. El producto escalar de un vector por sí mismo es A ■ A = \A\ 2 . 

4. El producto escalar es conmutativo: A ■ B = B ■ A. 

Si descomponemos los vectores en coordenadas cartesianas: 
A = A x % + Ayj + A z k, 

B = B X Í+B y j + B z k, 

entonces, por las características de los vectores unitarios, tenemos que: 

Á-B = A X B X + A y B y + A Z B Z 

2.2.3 Producto vectorial 

Es otro tipo de producto entre vectores. Mientras que en el producto escalar se 
obtiene como resultado un escalar, en el producto vectorial el resultado es un 
vector. La notación que se usa es un aspa (x). Multiplicamos los vectores A y 
B que teníamos antes. 

ÁxB = {A x ,A y , A z ) x (B x , B y , B z ) = 

= (A y B z - A z B y )i + (A Z B X - A x B z )j + {A x B y - A y B x )k 
Para calcularlo podemos utilizar las reglas del cálculo de determinantes: 

Ax B 

Entonces : 

• El resultado es perpendicular a A y a B. 

• Si A y B son paralelos, el producto vectorial es cero. 

• El módulo del producto es \A x B\ = |A||_B| siné>, con 9 el ángulo entre los 
vectores A y B. 

• Ax B = -Bx A. 



i 


j 


k 


1, 


Ay 


A 


% 


By 


B, 



2.3 Las leyes de Newton 



Antes de empezar definiremos lo que es un sistema inercial: 

Sistema Inercial: Si sobre un objeto no actúa fuerza alguna, cualquier 
sistema de referencia respecto al cual el objeto se mueve sin aceleración es un 
Sistema Inercial. Las leyes de Newton se refieren a sistemas de referencia iner- 
ciales. 

1. Primera ley: Toda partícula libre se mueve con v = cte. 

• Una partícula libre es aquella no sometida a ninguna fuerza. 

• Esta ley se suele denominar Ley de Inercia. 
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2. Segunda ley: Si una partícula de masa m se mueve con aceleración a, es 
porque sobre ella actúa una fuerza tal que: 



F = ma 



• La fuerza la medimos en Newton (N) en el Sistema Internacional. 

N = Kg • m/s 2 

• Esta ley supone una definición de fuerza. 

• Principio de superposición: Si sobre una partícula con masa m actúan 
varias fuerzas i*\ , F 2 , ■ ■ ■ , F n , la aceleración resultante con que se 
mueve es a = di + a 2 + ■ ■ ■ + a n , donde a¿ es la aceleración que 
m tendría si solo actuase la fuerza í¿. La fuerza neta resultante es 
la suma vectorial de cada una de las fuerzas: 

n 

F = ma = m(ai + a 2 -\ h a n ) = Fi + F 2 H h F n = J^ F, L 

3. Tercera ley: Si un cuerpo A ejerce una fuerza Fab sobre un cuerpo B, 
entonces B ejerce sobre A una fuerza Fba de módulo y dirección iguales 
a la anterior pero de sentido contrario. 

Fab - -F BA 

• Se la denomina Ley de Acción- Reacción. 

Existen 4 interacciones fundamentales, o formas básicas de interacción entre 
dos partículas 

• Gravitatoria. Debida a las masas de las partículas. 

• Electromagnética. Debida a las cargas eléctricas. 

• Nuclear fuerte. Cohesión de los núcleos atómicos... 

• Nuclear débil. Estabilidad nuclear, radiactividad... 

Toda interacción real se reduce a combinaciones de éstas. Sin embargo, en 
ocasiones es muy complicado quedarse con estas cuatro leyes y se recurre a 
descripciones fcnomcnológicas: 

• Gravedad en la superficie terrestre: 

F = mg 



• Fuerza normal: 

• Ley de Hooke (en ID): 



N 



F = -kx 

• Tensiones transmitidas por cuerdas: 

F = f 

• Rozamientos estático y cinético: 

\P\=He\Ñ\, \P\=»c\tí\ 
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2.4 Sistemas de referencia Inerciales y No Iner- 
ciales 

Supongamos que tenemos dos sistemas de referencia, uno de ellos inercial con 
respecto a la posición de la partícula P (sistema O) y otro no inercial, (sistema 
O'), que se encuentra en la posición R(t) con respecto a O, y se mueve con 
velocidad V(t) y aceleración A(t). 




La posición de la partícula P es, para O, r(t), mientras que para O', la 
posición de P es r'(í): 

O : f(t) = x(t)í+ y(t)j + z{t)k 

dx^ dy - dz -> 
v{t) = ——i -\ — — j H — — k 
v ; dt dt J dt 

d 2 x^ d 2 y^ d 2 z j> 



O 1 : r'(t) = x'{t)i' + y'(t)f + z'(t)k' 
, , . dx' -> dy' -> d'z ,-» 

., s d 2 x' - d 2 y - d 2 z - 

¿Cómo se relacionan r y r'l De la figura anterior deducimos que: 

f=R + P 
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Para ver cómo se relacionan las respectivas velocidades que mide cada ob- 
servador, el observador (O) deriva la expresión anterior: 

dr dR d ( , - , t¡ , ,~,\ 
¿T = — = — + — xV + j/Y + z'fc' = y + w' 
ai dt dt \ ) 

Al realizar esta operación hay que tener en cuenta que los vectores unitarios 
i' , j', fe' son fijos para ambos observadores (y por tanto sus derivadas son cero), 
porque aunque (O') se mueve con respecto a (O), sus unitarios siempre tienen un 
mismo módulo y señalan en las mismas direcciones. Otra cosa sería si hubiese 
una rotación de un observador con respecto del otro: entonces los unitarios 
de (O') sí que cambiarían de orientación para (O) y sus derivadas serían más 
complicadas. 

De nuevo, (O) deriva para relacionar las aceleraciones que cada uno mide: 

_. d 2 r d 2 R d 2 



dt 2 dt 2 



+ -rp (x'V + y'f + z'k'\ = Á 



Es decir, al pasar de un sistema de referencia (O) a otro (O'), y si el 
movimiento de uno con respecto del otro es de traslación (y no rotación) tal 
y como se ha descrito en la figura anterior, el cambio en los vectores r, v y a es 
tal que: 



f=R 



v = V + v' 



a = A + a> 



Esto se ha de tener en cuenta a la hora de aplicar las leyes de Newton. 



2.5 Dinámica de una partícula 

El estado dinámico de una partícula viene dado por una serie de datos: masa, 
posición y velocidad en cada instante de tiempo respecto a un O. 

Si conocemos la masa de una partícula, y esta es constante, el estado dinámico 
de una partícula se determina mediante la 2 a Ley de Newton, que nos permite, a 
cada instante de tiempo, conocer como cambian r(t),v(t) — > f(t + dt), v(t + dt). 

f(t + dt) = r(t) + df= f(t) + v(t)dt 

F(i) 

v(t + dt) = v(t) +dv = v(t) + a(t)dt = v(t) + —^dt 

m 

Así que si conocemos F(t), utilizando la 2 a Ley y lo que sabemos de cinemática, 
podemos conocer punto a punto el estado dinámico. Pero lo que hemos hecho 
ha sido obtener una ley diferencial, que nos permite pasar de t a t + dt. Si lo que 
queremos es tener una expresión analítica de f(t) y v(t), entonces necesitamos 
integrar las ecuaciones anteriores utilizando las condiciones iniciales fo y ^o : 

v(t) = v Q + — / F(t)dt 

m Jt 
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r(t) =f + v(t)dt 

Jt a 

Como este cálculo es en muchas ocasiones inviablc analíticamente, es útil in- 
troducir otras magnitudes dinámicas que nos permitan obtener información del 
sistema y simplificar los cálculos. 



2.5.1 Momento lineal 

Se define como: 

Se mide en Kg ■ m/s. 



p = mv 



Relación con la 2 a Ley : En casos de masa constante es muy útil, puesto 
que se relaciona fácilmente con la fuerza total Fj- que actúa sobre un sistema: 

=• _ dv d dp 

Ft = ma — m— = — (mv) = — 

dt dt di 

Relación con la I a Ley : Cuando no hay interacción ninguna, Ft = 0, 
gracias a la ecuación anterior deducimos que: 



F T = = 



dp 
dt 



p = cte 



2.5.2 Momento angular 

El momento angular de una partícula J, definido respecto a un origen de coor- 
denadas O, es: 



Jo = rx p 



donde O es el origen desde donde se mide r. Al cambiar de origen, cambia el 
resultado: desde O mido Jo, y desde O' mido Jo 1 , pero Jo ^ Jo' ■ 




El momento angular es una magnitud útil porque en algunos casos se con- 
serva. Vamos a ver cómo es su derivada: 



dJp 
dt 



dr dp 

Tt xp + rx Tt 



J = rxp 

Calculamos por separado los dos sumandos que nos han salido: 

Jo 



^ x p = v x mv = 0) dJ c 



rx a £=rxF 

dt 



di 



= rx F. 
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Entonces, si encontramos un sistema de referencia O en el que el vector posición 
r sea paralelo al vector fuerza F, entonces ^JP- — 0, por las propiedades del 
producto vectorial, y se conservará el momento angular. 
Se mide en Kg ■ m 2 / s. 

2.5.3 Momento de una fuerza 

El momento de una fuerza No, definido respecto a un origen O desde el que 
medimos la posición R del punto donde se aplica la fuerza F, se define: 

Ñ n =rxF 



El momento de una fuerza nos da cuenta de cómo de eficaz es la fuerza F 
a la hora de rotar un cuerpo con respecto a un eje. El ejemplo tradicional es 
el de una puerta. Aplicando la misma fuerza justo al lado de las bisagras no 
obtenemos el mismo resultado que si lo hacemos en el extremo opuesto de la 
puerta, junto al pomo. 

El momento de una fuerza es importante en estática. Para que un sistema 
esté en equilibrio es necesario que se cumplan dos condiciones, que el sumatorio 
(o la resultante) de todas las fuerzas que actúan sobre el sistema sea cero, y que 
el sumatorio (o la resultante) de sus momentos también sea cero. 

2.5.4 Impulso de una fuerza 

El impulso F de una fuerza F que actúa entre los instantes de tiempo ti y 
Í2 > ti, se define: 



I = i F(t)dt 



Como se vio con anterioridad, F = 3? => Fdt = dp 

r«2 



/•Z2 

1= / dp = p(t 2 )-p(t 1 ) = Ap 
Jti 



es decir, / nos da información que depende de los estados inicial (en t = ti) y 
final (en í = f 2 ). 

Es útil porque nos permite conocer cuál ha sido la fuerza media que se ha 
estado aplicando entre ti y ¿2, 

También, en caso de ser nulo, nos permite afirmar que p(t 2 ) = p(t\). 

2.5.5 Trabajo 

El trabajo diferencial realizado por una fuerda F que actúa sobre una partícula 
que se desplaza una distancia dr se define así: 

dW = F -dr = F dr eos 9 
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f(t + dt) 



En una trayectoria que va desde una posición inicial f\ hasta una posición 
final r*2, integramos la expresión anterior para conocer el trabajo total realizado: 



W = F-dr 

Jri 



Si la fuerza es constante la puedo sacar de la integral: 



W = F ■ dr = F- (Af) 



Las unidades de medida del trabajo en SI son los Julios: 1J = 1N x m. Un 
Julio equivale a un Newton de fuerza por metro. Por la definición de N: J — 
Kg m 2 /s 2 . 



2.5.6 Potencia 

La potencia P se define como el ritmo con el que una fuerza F realiza un trabajo, 
es decir, como trabajo W por unidad de tiempo: 



P = 



dW 
~dT 



como dW — F ■ dr = F ■ vdt, sustituimos en la expresión anterior: 



P = F -v 



Esta es la potencia instantánea. Entre los instantes de tiempo t\ y ¿2, la potencia 
media (P) se define como: 



(P) = 



1 

Á7 



P(t)dt = 



1 

Af. 



t2 W 

dW = — 

Ai 



La potencia se mide en vatios: 1W = U/s. Un vatio es un trabajo de un 
Julio realizado en un segundo. 
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2.5.7 Energía cinética 

La energía cinética de una partícula que tiene masa m y se mueve con velocidad 
v, se define: 



n 1 2 p 2 

E r . = -mv = 

2 2m 



Las unidades son: Kg m / s = J los Julios. Misma unidad que el trabajo. 
Es también habitual utilizar la letra T. 

Teorema Trabajo-Energía Cinética 

Recordemos la expresión del trabajo realizado por una fuerza que actúa sobre 
una partícula entre sus posiciones ri yr^: 

rr 2 

W= F-df 

Jri 

Como df = vdt, y F — ma = rn%, escribimos: 

ñ .. _ ñ dv 

W = F ■ dr = m— • vdt 

J f\ J r\ 

Ahora nos aprovecharemos de la relación diferencial: 

d(v ■ v) dv 

; = ¿v ■ — 

dt dt 

y por lo tanto el trabajo nos queda: 

r 2 m 1 

W = J^ -d(v 2 ) = -m{v 2 2 - v\) = E c2 - E cl = AE C 



W = AE r 



El cambio de energía cinética en una partícula durante su trayectoria es igual al 
trabajo que las fuerzas que actúan sobre ella han realizado a lo largo de dicha 
trayectoria. 

Las fuerzas no siempre realizan trabajo. El ejemplo típico es el de movimiento 
circular uniforme, con aceleración centrípeta, dirigida hacia el centro, y por lo 
tanto la fuerza también. El diferencial desplazamiento es siempre perpendicular 
a la fuerza centrípeta, por tanto: F _L df =>■ W = 0. 

2.5.8 Fuerzas Conservativas 

Como hemos visto antes, el trabajo total realizado por una fuerza entre dos 
puntos r*i y fi es una integral. Y por lo tanto, en general, depende de cómo se 
haga ese camino desde la posición f\ hasta la r^. Pero algunas fuerzas, llamadas 
conservativas, se caracterizan por realizar un trabajo que solo depende de las 
posiciones inicial y final, y no de lo que pase durante el camino. Dicho de 
otro modo, son fuerzas que, si actúan durante una trayectoria cerrada, donde el 
punto final y el inicial coinciden, realizan un trabajo total cero: 

F r -dr = 
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para todo camino cerrado C . 

Como se ha dicho antes, el trabajo en una trayectoria cualquiera, entre r\ y 
r 2 es: 

rr 2 

W(n ->r 2 )= F-dr 

J r i 

y si la fuerza es conservativa, esta cantidad solo depende de r\ y r 2 =>■ si 
escogemos un punto de referencia r IO ¡, el trabajo que una fuerza realiza para 
trasladar una partícula desde un r general hasta la referencia escogida es una 
cantidad que solo depende de r. Esto podemos decirlo porque trabajamos con 
fuerzas conservativas. 

2.5.9 Energía Potencial 

Entonces, como el trabajo para traer nuestra partícula hasta la referencia f re f 
solo depende de r, podemos definir la energía potencial, una función escalar que 
dependa de r como esa cantidad de trabajo realizado: 



E p {r) = W(f^ f re f) = / F(r) ■ dr = F{r) ■ dr 

Jr Fr loí 

Normalmente lo más cómodo es asignarle a la r rc f una energía potencial cero. 
Entonces, la energía potencial en r, es la cantidad de trabajo que haría falta 
para llevar una partícula desde r hasta r re f, o dicho de otro modo, la cantidad 
de trabajo, con signo cambiado, que se habría tenido que realizar para llevar la 
partícula desde r re f hasta r. 

Al llevar una partícula de una posición f\ hasta r 2 , y como estamos con 
fuerzas conservativas, se hace el mismo trabajo independientemente de la ruta 
seguida: El trabajo total será el mismo si seguimos la ruta de puntos o si vamos 
por la línea discontinua, es decir, por el tramo T\ y luego por el tramo T 2 . 




r 2 ) = W(n -> f ref ) + W(r rcí -^ r 2 ) 
r rcf ) = E p (fi) 

W(r rcf -». r 2 ) = -E p {r 2 ) 

Entonces podemos seguir la ruta propuesta en el esquema superior: para ir 
de r\ a f 2 podemos pasar por f re f. El trabajo realizado es 

W(n -4 r 2 ) = W(n -4 r ref ) + W(f Tei -4 f 2 ) 

De estos dos sumandos, el primero es, por definición, E p {f\), y el segundo es 
—E p (f 2 ). Entonces, la cantidad de trabajo que se ha realizado para ir desde f\ 
hasta r 2 es: 

W(n -> f 2 ) = Epin) - E p (f 2 ) = -AE P . 
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Es decir, tal y como hemos definido la E p , podemos afirmar que, con fuerzas 
conservativas, el cambio de energía potencial al desplazar una partícula de una 
posición a otra es igual al trabajo realizado por dicha fuerza con signo negativo. 

W = -AE P 



Para la energía potencial a veces se utiliza la notación U. 
2.5.10 Energía Mecánica 

Recordatorio: 



Con F conservativa: W = —AE p 



Con cualquier F: W = AE, 

vi 

A(E C + E P )=0 

? p se mantiene c 
tivas. Esta cantidad se denomina Energía Mecánica. 



Con F conservativa: AE C = —AE p 



Es decir, la cantidad E c + E p se mantiene constante con Fuerzas conservá- 



is — E c + E p 



Utilidades de la Energía Mecánica 

La energía proporciona información importante: Si conocemos el valor de la 
energía mecánica en un punto concreto fg , y estamos en un sistema conservativo, 
ya sabemos que esa cantidad se va a conservar: la energía mecánica total no 
cambia, y vale lo mismo en cada posición f: E(fo) = E(r). 

Si además conocemos la expresión de la energía potencial E p (r), y las condi- 
ciones iniciales, podemos determinar el movimiento: 

1 1 / dr\ 2 í dr\ 2 9 

E^-mv^E^E^-m^-) =E ^ (-) =-(£-2^) 

df r {t) df 

dt = — . => t — í 



^(E-E p (r)) ' Jro ^ {E -E p {r)) 

de esta última expresión podemos despejar f(t). 

2.6 Ejemplos de Fuerzas Conservativas 

2.6.1 El caso Gravitatorio 

Una partícula de masa mi situada en la posición r\ ejerce sobre una masa mi 
ubicada en la posición r 2 , una fuerza : 

^ Gm\m 2 _ 
F\2 = ^ r i2 

Donde el vector r*i2 es el que va desde la partícula 1 hasta la partícula 2: 

r\2 = r 2 - n 

por otro lado G es la constante de gravitación universal, y vale G = 6, 67 ■ 
1CT 11 Nm 2 /Kg 2 . 
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2.6.2 Energía potencial gravitatoria 

Antes de establecer una expresión para la energía potencial gravitatoria debe- 
mos comprobar si la fuerza de la gravedad es conservativa. Es decir, debemos 
comprobar si el trabajo para ir desde una posición r¿ inicial a una final ff 
depende de la trayectoria seguida. 



W{n -+ f f ) 



F\2 ■ df = —Gm\mi 



f 12 ■ df 






fu(í + dt) 



del esquema superior podemos deducir que: 

f 12 -df= |ri 2 ||dr|cos0 = |fi2|d|fi 2 | 



nidn 



porque |<ir|cos# es la proyección del vector dr sobre ri2, y coincide con la 
variación del módulo de éste último. Obsérvese que \fi2(t + dt)\ — \fi2(t)\+dri2- 
entonces, el trabajo realizado lo escribimos como: 



W(fi — > ff) = —Gmim,2 



dri2 
ri2 2 



= -\-Gm\rri2 



1 
ri2 



= Gm\m2 



1 



1 



es decir, el trabajo que realiza la fuerza de la gravedad cuando la partícula 
va desde r¿ hasta ff no depende del camino seguido sino de dichas posiciones 
inicial y final. Por lo tanto es una fuerza conservativa y podemos definir la 
energía potencial gravitatoria de un punto f general a partir de un f rc ¡ que 
escojamos: 



E P {r) = - / F12 ■ df= -Gmim,2 



1 



1 



/ r rof \ ' r rof 

Por comodidad aquí conviene tomar la referencia r rc ¡ 

mim.2 



= E p (r) 



00, asi: 



E p (r) = -G- 



Los cuerpos de masa M y simetría esférica generan, fuera de la esfera de 
radio R, un campo gravitatorio igual al que existiría si en su centro hubiese una 
partícula puntual de la misma masa. 

Cerca de la superficie, es decir, a altura h (distancia del centro r — R + h) 
pequeña (r = R + h ~ R), la atracción sobre una masa m es: 



F = -G^u r 



mg 



Fundamentos de Física I 



Guillermo Sánchez 



Así definimos la g, aceleración del campo gravitatorio, aproximadamente cons- 
tante si nos movemos en distancias a la superficie que sean muy pequeñas com- 
parada con el radio terrestre. 

La fuerza del peso es F = mg = —mgk. Entonces, para estos casos en 
que tenemos un campo uniforme podemos volver a calcular la energía potencial 
gravitatoria: 



E P (f) 



F-dr 



-mg 



(—k)-dr = mg / dz — mg(z — z Te {) — E p {z) 



Si escribimos la altura respecto a una referencia como h = z — z rc ¡ tenemos la 
típica expresión: 



E p = mgh 



Recordar que esta expresión vale solo para campos gravitatorios creados por 
una masa sobre volumen esférico en las proximidades de su superficie. 

2.6.3 El caso elástico 

En un resorte ideal la fuerza recuperadora es proporcional a la deformación. 
Supongamos que el muelle, cuando no está deformado, tiene una longitud na- 
tural Iq y mediante la acción de una fuerza externa F cxt lo deformamos una 
distancia x hasta tener una longitud l = Iq + x. Entonces, según la Ley de 
Hooke, aparece una fuerza F r en el muelle tal que: 



F« 



-F r = k(l - l )i 



-kx 



lo 





l 

donde k es una constante. 

Conviene recordar que la fuerza recuperadora siempre es en el sentido opuesto 
a la deformación. Si la deformación es positiva, la fuerza es negativa y viceversa. 

Vamos a comprobar a continuación si la fuerza recuperadora es conservativa. 

El trabajo realizado al arrastrar una partícula sujeta al extremo de un muelle 
desde x\ hasta X2 es: 



W(xx -> x 2 ) = 



F r ■ dx = 



(—kxi) ■ (dxi) = — k 



X2 -| 

XQiX — K\ X'y .Ti I 



Xl 
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El resultado como vemos es una función que depende de los extremos inicial y 
final del trayecto seguido, pero no del camino. Por tanto es una fuerza conser- 
vativa y podemos definir una energía potencial elástica. 



E p {x) 



/ F r dx= -(x 2 -x 2 e{ ) 

Jx rcf z 



Se escoge como referencia el origen de coordenadas, x Te { = 0, de manera que 
E p (x) = \kx 2 , con x = l — lo- 

Entonces, podemos afirmar que si sobre una partícula de masa m actúa 
únicamente una fuerza elástica de Hookc, su energía mecánica se conserva: 

E = -mv H — kx = cte 
2 2 

2.7 Dinámica de Sistemas de Partículas 

Llegados a este punto pasaremos a considerar sistemas de partículas: conjuntos 
de partículas en el estudio, que interaccionan entre sí y con otros sistemas. 

Antes de nada definimos el centro de masas. Si hay n partículas, cada una 
de Has con una masa m¿ y posición r¿, entonces, la masa total es: 

M = ^2 m i 

i— 1,71 

y el centro de masas R: 

MR= ^2 m i^ 

i—l.n 

Más tarde nos será útil para entender las interacciones presentes en el sis- 
tema. Veamos como son dichas interacciones. 

íoxt 



O 

En el esquema hay dibujadas n partículas. La ¿-csima, o partícula i interac- 
ciona con todas y cada una de las demás partículas que conforman el sistema 
(Fuerzas Internas) y por otro lado está sometida a una Fuerza Externa o resul- 
tante de las fuerzas externas -F e xt ¡. 
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Obsérvese que las fuerzas internas las hemos denotado de la siguiente ma- 
nera: La fuerza que actúa sobre i como consecuencia de su interacción con la 
partícula 2 es F i2 etc.. 

Entonces la segunda ley de Newton, para esta partícula la escribiríamos: 

j2- &* 



— 2_^i *3 """" c: 



midi = w¿ — 

.7 = 1,71 

Ahora aplicamos la segunda ley al sistema de partículas en su conjunto, es 
decir, sumamos la expresión anterior para las n partículas. Pero antes podemos 
tener en cuenta un aspecto de las fuerzas internas, y es que, debido a la ley de 
acción reacción, por cada fuerza interna abrá una igual y de sentido contrario: 
Para cada pareja de partículas i y j se anulan las fuerzas internas entre ellas: 

hitó ( _^ 

Y miSi = Y [ Y ñj ] + Y Fext .¡ 

i—l,n i \ j I i—l,n 




Entonces, la segunda ley de Newton la podemos reescribir así: 

y J mÁ — 2_^ ^ext.i — ^cxt,T 
¿=l,n i=l,n 

Y m í«¿ ^Y mi ^ff = ^ {Y mi ^ i ) = ^oxt,T 

La fuerza externa resultante que actúa sobre el sistema, F e *t,T, es la derivada de 
la suma de los momentos lineales individuales de cada partícula. A partir de esta 
expresión cobra sentido la definición de momento lineal del sistema: P — ^Pi- 
De modo que obtenemos una ecuación análoga a la que ya conocíamos para una 
partícula, que ahora es válida en sistemas de más de una partícula. 



F - dP 
dt 

Al igual que en una partícula escribimos p = mv, en un sistema podemos 
escribir P = MV. La masa M es la masa total del sistema, pero, ¿Cuál es la 
velocidad VI . 

s% _ P _ Y. m iVi _ d (Y, m i^i 



M J2 m i dt \ J2 m 

es decir, V es la velocidad con que se mueve el centro de masas R que habíamos 
definido al comienzo de este apartado. 

Recapitulando, tenemos que, en sistemas de varias partículas: 

P = MV 

F ext , T = § = M§ = MA 
dt dt 
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son expresiones análogas a las de una partícula, o más bien generalizaciones. 

De estas ecuaciones se deduce que el movimiento del centro de masas es 
independiente de cómo sean las fuerzas interiores, y tampoco depende de la 
fuerza exterior aplicada a cada partícula del sistema: solo depende de la fuerza 
exterior resultante. 

Además, al igual que sucedía con una partícula, en caso de fuerza externa 
resultante nula, se conserva el momento lineal: 

dP _ - 

—¡T — -fcxt.T 

dt 
si ^oxt.T = =*> P = cíe 

2.8 Colisiones entre partículas 

2.8.1 Partículas indeformables y libres 

Estudiamos el caso de colisiones entre esferas indeformables y aisladas, no 
sometidas a fuerza alguna: 

mi, Vi m 2 ,v 2 



mi, Mi m 2 ,u 2 



Como hemos dicho que están aisladas, no existen fuerzas externas, y por lo 
tanto el momento lineal del sistema se conserva: 
Antes del choque, tenemos: P¿ = miVi + m 2 v 2 . 
Después de la colisión: Pf — miUi + m 2 u 2 . 

m\V\ + m 2 v 2 = miüi + m 2 u 2 (2.1) 

Además suponemos que la colisión es frontal, todas las velocidades, antes y 
después, apuntan en la misma dirección (aunque los sentidos puedan ser opues- 
tos), la que marca el vector unitario i. 

Como se ha dicho antes, las esferas no se deforman, así que no se realiza 
ningún trabajo interno (de deformación). Esto unido a la ausencia de fuerzas 
externas nos permite afirmar que no hay trabajo ninguno durante el proceso: 
W = => E c = cte. 

Antes del choque la energía cinética es E ci = \m\v\ + \vn 2 v\. 

Y después: E c ¡ = \m\%i\ + \m 2 u 2 . 

Entonces: 

-m\v\ + -m 2 v\ = -miWi + -m 2 u 2 2 (2.2) 

A partir de las ecuaciones 2.1 y 2.2, y suponiendo conocidas las masas y 
velocidades iniciales, deducimos las velocidades finales. 

(mi - m 2 )ui + 2m 2 v 2 

ui = ■ 

mi + m 2 
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(m 2 — mi)v2 + 1rri\V\ 

u 2 = . 

mi + m 2 

Cuando se conservan tanto el momento lineal como la energía cinética se 
dice que se trata de una colisión elástica. 

La velocidad relativa entre las partículas se conserva: u\ — u 2 = v 2 — Vi ■ 

2.8.2 Colisión no rígida 

Si la colisión es no rígida, con deformación, se puede realizar un trabajo interno 
W¡nt durante la misma, en el caso de que la fuerza interna sea no conservativa. 
En ese caso la energía cinética no se conserva. Se emplea entonces el llamado 
coeficiente de restitución / que nos señana cómo ha sido la variación de la 

velocidad relativa: 

_ u 2 - U\ 

Vi - v 2 

Si conocemos / y las velocidades iniciales v\ y v 2 , podemos despejar las 
finales: 

m\V\ + m 2 v 2 — m\U\ + m 2 u 2 \ 

«2-Wl = /(í>2-Vl) J 

2.8.3 Colisión no aislada 

Ahora lo que sucede es que el momento no se conerva, porque existen fuerzas 
externas. 

F ext + => P ¿ cte 

Como sabemos, la variación del momento lineal la podemos conocer a través 
de la variable Impulso: 

/■Ai 
/ext = / F cxt dt = AP 

Jo 

Si la colisión es percusiva, es decir, si Ai ~ 0, la integral de una fuerza finita es 
aproximadamente cero, así que sigue siendo buena la aproximación de P¿ ~ Pf, 
tomando los instantes inmediatamente anterior y posterior al choque. 
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Chapter 3 

Oscilaciones 

3.1 Introducción 

Entendemos por oscilación el movimiento que surge tras la perturbación de un 
sistema desde una posición de equilibrio estable. La oscilación es un "vaivén" 
en torno a una posición de equilibrio. 

El modelo más sencillo es el de la Oscilación Armónica, que ya conocemos. 
En el caso de que la posición de equilibrio sea x CCi — 0, el sistema viene descrito 
por la ecuación diferencial característica: 



x 



,2 



üj ¿ x = 



dt 2 

3.2 Cinemática de la oscilación armónica 

La solución a esta ecuación es una función sinusoidal: 

x(t) = A sin(cjí + a) = A cos(u>t + ¡3) , ¡3 = a — — 

La amplitud es A, la fase inicial es a (o ¡3). La velocidad y aceleración se 
determinan por simple derivación: 

v(t) = — = Alo cosíwí + a) 
dt 

. , (X Jb\Jj \ , O . / \ Q 

a(t) = — -t. — = —Au> srníwí + a) = —lo x 
dt z 

Estas funciones se caracterizan por su periodicidad. Pasado un periodo P: 
x(t) =x(t + P)\ ^ p= 2tt 



v(t) =v(t + P) 



Lü 



La frecuencia es v = 1/P, y la frecuencia angular es la ui de la ecuación diferen- 
cial característica u = 2ir/P = 2ttv. 

En caso de que la posición de equilibrio sea otra: a; cq ^ 0, cambia algo la 
ecuación característica, porque ahora la posición es: 

x(t) = x cq + A sin(wí + a) =^ 
39 
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(Ai JL\L j a "2 • / \ 2 / \ 

— — - — = — Alo sm(wí + a) = —lo (x — x cq ) 

Es decir, en general, con x cq ^ 0, el movimiento oscilatorio armónico se rige 
por la ecuación diferencial característica: 

2/ \ / c\ -t \ 

— + lo (x-x cq ) (3.1) 



Conociendo las condiciones iniciales podemos deducir amplitud y fase ini- 
cial. Es decir, conociendo posición y velocidad en el instante t — 0, (xo y vq 
respectivamente), podemos averiguar Aya: 

xq = x(t = 0) = x cq + A sin al ( xq — x cq = Asina 

Vo = v(t = 0) = Alj eos a J l 77 ~ Acosa 

Dividiendo estas dos últimas ecuaciones: 

u>(x - x cq ) 
— = tan a 

Mientras que elevando al cuadrado esas mismas ecuaciones y sumando: 

En resumen, conociendo xq y vq tenemos que: 

A = J(x a - x cq ) 2 + í — 

Lj(x - X cq ) 



a = arctan 

3.3 Dinámica de la oscilación armónica 

Ahora pasamos a la dinámica de este tipo de oscilaciones — > nos preocupamos 
de la fuerza que la produce, su energía... 

Utilizando la ecuación x(t) y la seguna ley de Newton calculamos la fuerza: 

x = x cq + Asm(Ljt + a)\ 2 . 2 

i, > => b = —muj A smlcjí + a) = —mu) (x — x ca ) 
F = ma = m^¡§ J 

Es decir, nos ha salido una fuerza que es proporcional a — (x — x cq ). Es 
decir, una fuerza de la misma naturaleza que la de los muelles. La constante de 
proporcionalidad mw 2 haría las veces de constante recuperadora del muelle k 

k = mu u> = \ikjm 

Cuando x = x cq , la fuerza es cero F(x eq ) = 0. En caso contrario es propor- 
cional al desplazamiento x — x cq dirigida hacia x cq . 

Como es del mismo tipo que la fuerza en muelles, es decir, sigue la ley de 
Hooke, sabemos que es una fuerza conservativa (como se vio en el tema 2). 
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Entonces, podemos definir una energía potencial E p . Para ello seleccionamos 
una referencia. En este caso lo más sencillo es tomar x rc f = x cq , para así tener 



KyX X Gq ) 



f X 1 

E p {x) = - -k(x - x c(1 )dx' = . . . - 

Jx Ia{ =x ac¡ 

Como vemos, mientras que la fuerza tenía dependencia lineal con el desplaza- 
miento, la energía potencial tiene una dependencia cuadrática. Esta es una de 
las características de la oscilación armónica. 
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La energía mecánica, como en todas las fuerzas conservativas, se conserva. 

E m = E c + E p 
E p (x) = ñk(x — x cq ) = — kA sin (utt + a) 

E c (x) = -mv 2 = -muí 1 A 2 cos 2 (u)t + a) 
Por tanto, si recordamos que k = muí 2 y sumamos las dos últimas expresiones: 



E m = E p (x) + E c (x) = hA 2 



La energía mecánica es constante. 



3.4 Péndulo simple 



Por péndulo simple entendemos aquél que tiene una masa puntual sostenida 
por una cuerda inextensible de masa despreciable, que no tiene pérdidas por 
rozamiento con el aire o con la pared de la que está suspendido. 
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A continuación tenemos un pequeño esquema con las fuerzas que actúan 
sobre la masa m cuando se desplaza un ángulo 9 con respecto a su posición de 
equilibrio. Actúan por un lado el peso y por otro la tensión del hilo. 




mg 
La fuerza resultante por tanto será la suma vectorial de las dos que actúan: 

F = T + mg 

Descomponemos las fuerzas en las direcciones radial y angular: 

F = mg sui9(—Üg) + (mg eos 9 — T)u r 

En la dirección radial la fuerza resultante es cero (notar que el cordel tiene 
longitud fija, l). En la dirección angular estamos fuera de equilibrio, con lo cual 
la suma de las fuerzas: 

F g = mag 



La aceleración angular o polar 



Fg = mag = mi 



es, para radio fijo: ag 
d 2 * 



, d 
dt 2 



. Por lo tanto, 



dt 2 



-mg sin 9 



despejando, obtenemos la siguiente ecuación del movimiento: 



dP 



9 ■ a 
— sint* 







En el caso de que las oscilaciones sean de pequeña amplitud, es decir, con 9 
bajo, podemos aproximar el seno por el ángulo. 

Si 9 — > =^ sin# ~ 9 (¡Ojo! trabajamos con radianes). 

Entonces, la ecuación para pequeños ángulos es la de una oscilación armónica: 



dt 2 



I e = ° 



El ángulo de equilibrio es obviamente 9 e 



0. 



Comparando esta ecuación que hemos obtenido con ecuación diferencial ca- 
racterística del movimicno oscilatorio armónico (ecuación (3.1)), tenemos que 
la frecuencia es ui 2 = g/l, mientras que el periodo es P = 2w^/l/g. Además, la 
evolución del ángulo 9 responde a una ecuación del tipo: 9(t) = 9osm(ujt + a), 
con amplitud 9q. 



- Fundamentos de Física I - Guillermo Sánchez 

3.5 Oscilaciones amortiguadas 

Hasta ahora hemos estudiado el oscilador armónico ideal, 



d 2 x , ... . 

m—rr = —kx =>■ x = AsmUjjQt + a) 
dt z 



y con frecuencia natural u>q = yjk/m. 

Sin embargo es difícil encontrar en la naturaleza situaciones ideales. Los 
osciladores más realistas son aquéllos en los que actúa también una fuerza disi- 
pativa, algún tipo de fricción que se opone al movimiento. El modelo más 
comunmente empleado para describir estas situaciones es el de una fuerza pro- 
porcional a la velocidad del oscilador, pero con sentido contrario a ésta: 



F = —cv 



con c =cte. Entonces a la ecuación característica le tenemos que añadir un 
término: 

d 2 x dx 

m—rr = —kx — c—- 

dt 2 dt 



como k/m = u!q, 



d 2 x c dx 2 



Definimos el coeficiente de amortiguamiento, 7 = ^~ La solución a este sis- 
tema, para amortiguamientos pequeños (podéis hacer el ejercicio) es la siguiente 
ecuación: 



j(í) = A exp 7t sin(w a í + a) 



donde la frecuencia de oscilación es 



^ = V w o - 7 2 



El resultado es una oscilación, con frecuencia menor que la frecuencia natural 
uj a < u>0) cuya amplitud está modulada por una envolvente decreciente. Es 
decir, la amplitud efectiva está amortiguada exponcncialmcnte A{t) — A exp -7 * 
(líneas de puntos en el gráfico). El resultado es válido sólo para amortiguamiento 
pequeño, es decir, para 7 < cjq- 
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10 



15 



20 



25 



A mayor 7, más rápido será el amortiguamiento. A la inversa del coeficiente 
de amortiguamiento se le llama tiempo de relajación r = I/7. 

Por otro lado el período de las oscilaciones es P a = 2ir/uj a . 

Las cantidades A y a, es decir, amplitud y fase inicial, se pueden obtener si 
se conocen posición y velocidad iniciales Xo, Vq. 



3.6 Oscilaciones forzadas y amortiguadas 

En este apartado damos una vuelta de tuerca más al oscilador armónico. Si, 
como sucede de forma natural, hay algún tipo de amortiguamiento, pero que- 
remos que la oscilación se mantenga, entonces necesitamos aplicar al oscilador 
una fuerza impulsora periódica: F — Fq cos(wí) 
Entonces, la ecuación dinámica queda: 



d*x 



-kx 



dx 



Fq COs(wí) 
F. externa impulsora 



dt 

Amortiguamiento 

dividiendo por m y sustituyendo: k/m = Wq, c/m = 2j 

d 2 x 2 dx F 

_ +Woa; + 27 _ = _ cos(u;í) 



(3.2) 



Es decir, tenemos una ecuación que es por un lado una mezcla de oscilación 
amortiguada (~ exp -7 *, a frecuencia ui a = y Wq — 7 2 , y por otro estamos 
forzando al sistema a oscilar con una frecuencia u, de la fuerza impulsora. A la 
frecuencia de oscilación natural la llamamos ujq. 

Cuando el sistema comienza a funcionar, pasa por una fase transitoria 
inicial. El oscilador tiene una frecuencia inicial uj a (la amortiguada). Pasado 
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un tiempo mucho mayor que el tiempo de amortiguamiento, t >> t = 1/7, 
la energía que pierde el sistema (por el amortiguamiento) y la que recibe por 
parte de la fuerza impulsora se equilibran. Entonces se alcanza la solución 
estacionaria. La solución del sistema es del tipo: 



x(i) = B cos(urt — tp) 



El sistema oscila con la frecuencia forzada lo. Veamos cuanto valen la ampli- 
tud B y la fase tp: En primer lugar calculamos las derivadas primera y segunda 
de x con respecto a t: 

dx 

= —Bujsm(üjt — tp) 



dt 

d^x 

~d¥ 



-Bui cos(wí — tp) 



Sustituimos las expresiones de x, dx/dt y d 2 x/dt 2 en la ecuación característica 

(3.2): 



d x 2 dx Fq 

_ +WoX + 27 _ = _ cosM) 



F 



—Bui cosÍLot — ip) + lüqB cosíujt — tp) — 27 — Bu sin(cjí — tp) = — cos(cjí) 

m 

Recordemos dos relaciones trigonométricas: 

cos(wí — tp) = cos(wí) eos tp + sin(u;í) sin tp 

sin(cjí — tp) = sin(cjí) eos tp — cos(ujt) sin tp 
Sustituyendo en la expresión anterior: 

B (u! Q — cu ) (cos(wí) eos tp + sin(wí) smtp)—2B^ui (sin(wí) eos tp — cosíuit) sin tp) = — cos(wí) 

m 

Agrupamos sacando factor común cos(u;í) por un lado y sin(wí) por el otro: 

F ~ 



cos(wí) 



B eos tp (cjq — lo 2 ) + 2i?7 sin tp — 



+sin(u;í) \B sin tp (ui — lo ) — 2B-fU! eos tp\ = 



Cuando el coseno sea cero, el seno no lo será, y viceversa. Para que la 
expresión previa se cumpla para todo í, es necesario que sean cero tanto el 
factor que multiplica a cos(wí) como el que multiplica a sin(cjí), es decir, es 
necesario que sea cero lo que va dentro de los corchetes: 

Fn 

cos(cjí) — > B eos tp (uin — lo 2 ) + 2B-jlo sin tp =0 =^> 

m 



B 



F /m 



eos tp(tVQ — lo 2 ) + 27a; sin tp 
sin(wí) —¥ B sin tp (wq — u 2 ) — 2Bjlo eos tp = 

2"/LO 

tan tp = — 2 j 



(3.3) 



(3.4) 
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Ya tenemos el ángulo tp. Ahora la amplitud, pero antes recordamos otras 
dos relaciones trigonométricas: 



sinos = 



coso? = 



taníp 



\J\ + tan 2 os 

1 

\J\ + tan 2 p 



como hemos calculado la tangente, sustituimos y así tendremos el coseno y el 
seno despejados: 

2jlo 



sinos 



v /(w 2 -cj 2 ) 2 + 47W 2 



eos os = — = 

V(^-W 2 ) 2 +4 7W 2 

y por lo tanto, introduciendo estos valores en la ecuación (3.3), tenemos que: 

F /m 



B = 



v /(^-^ 2 ) 2 +4 7 2 w 2 



en resumen, tenemos un oscilador que, en estado estacionario, tiene frecuen- 
cia lo (que es la frecuencia a la que le forzamos a trabajar mediante la aplicación 
de una fuerza externa), fase ip y amplitud B: 



(3.5) 



B = 




F /m 


VK 2 - 


- lo 2 ) 2 +47 2 w 2 












taños = 


27W 

LO 2 — LO 2 





(3.6) 



Observar que ahora la amplitud y fase alcanzadas en el estado estacionario no 
dependen de las condiciones iniciales, sino que dependen de las características 
propias de la partícula oscilante (su masa m), la componente amortiguadora 
(wo, 7) y la fuerza externa impulsora (Fq, lo). 

Es interesante estudiar como cambian os y sobre todo B con la frecuencia 
forzada lo. 



B(lo) 



F /m 



a/(w 2 -cj 2 ) 2 + 47 2 w 2 

• si lo — 0, entonces B = B a = °' 2 m = F /k. 

• si lo — )• ex), entonces B — >• 0. No hay oscilación. 

• existe un valor de lo en el que la amplitud es máxima. 

— Máximo local en lo tal que dB/dLo = 
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,IB F " l [(lo 2 lo 2 ) 2 + 4jW] - 3/2 [2(oj 2 u 2 )(-2oj) + 87^] 



duj m 2 

dB 

duj 



U) = U¡ r = y Wq — 27 



es la llamada frecuencia de resonancia, aquella para la que la amplitud con- 
seguida es máxima. Si estimulamos el sistema con una fuerza que tenga esa 
frecuencia, la amplitud obtenida será: 



B„ 



B(iü r ) 



F /m 



F /m 



^2-^)2+4^2 2 7 ^4^f 



si recordamos, Fo/m = Bqüj 2 ,. Por otro lado uj a = yWg — 7 2 • Entonces la 
expresión se puede transformar para que quede así: 



Bmax — -Sn 



27^a 



Y gráficamente queda como se muestra en la figura: 




co,. 



co 



3.6.1 Resonancia 

Cuando se entra en resonancia se producen oscilaciones de gran amplitud, mucho 
mayores que la amplitud asociada a la fuerza externa, Fq — kB(¡. 
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Obsérvese que las amplitudes máximas se obtienen cuando se trabaja con 
coeficientes de amortiguamiento muy bajos en comparación con la frecuencia 
natural iúq. 

Wn 



S ° / n > B o «*• 7 < ^o 

2 7\/wo - 7 2 



El hecho de escoger valores pequeños de 7 tiene algunas implicaciones. En 
concreto, la frecuencia de resonancia se acerca mucho a la natural: 



us r = y w 2 - 27 2 ~ u} 
lo mismo sucede con la frecuencia de oscilación amortiguada: 



w o - 7 2 - w o 



entonces, la amplitud máxima queda así: 



Bq- ~ B — = B Q 



donde hemos definido el factor de calidad Q — lüq/2j. Cuando se entra en 
resonancia con 7 < wo, se tienen valores Q > 1 y por lo tanto -B max S* Bq. 
Además la frecuencia a la que oscila el sistema es muy próxima a la frecuencia 
de oscilación natural del armónico: uj r ~ lüq. 

La intensidad de la resonancia se asocia por tanto con la altura del pico. La 
curva de resonancia también queda caracterizada por la anchura del pico. Se 
suele tomar como referencia para calcular dicha anchura a las frecuencias a las 
que la altura es £? max /v2, frecuencias angulares que denotamos como w_ y u>+. 
Las frecuencias uj± son las soluciones de la ecuación de segundo grado que se 
obtiene al igualar la ecuación 3.5 a la amplitud deseada, en este caso -B max /v2, 
es decir: 



B(«) 



Fo/m 



V / K-cj 2 ) 2 + 47 2 w 2 V2 
Se obtiene una ecuación cuadrática cuya solución es, para 7 -C Wo, 



u(B = -B max /V2) = lüq ± 7 <! _ ° i- =^ Acj = w + - w_ = 27 

w + — cj +7 . 
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20 



15 



03 10 



Q 3 > Q 2 > Q! 






Qi - 

Q 2 

Q 3 




Bmax'« 




Anchura: co + -ü)_ - 2y l 


\ Q»l => (0 r = Cüg 


■\\ 




<B + 



10 



3.7 Oscilaciones Anarmónicas 



En las oscilaciones armónicas veíamos oscilaciones en torno a un punto donde 
la fuerza es cero (el punto de equilibrio x cq ), y una dinámica tal que la fuerza es 
proporcional al desplazamiento con respecto a dicho equilibrio: F oc — (x — x cq ) 
(ver gráfico derecho de la página 3). 

En un caso más general la fuerza no es lineal con el desplazamiento. Entonces 
no se obtiene una oscilación armónica en general, aunque es posible que, en 
torno a un punto de equilibrio, el sistema si que oscile. Un punto de equilibrio 
es cualquier punto x% donde F{xí) — 0. Pero para que haya oscilaciones dicho 
equilibrio debe ser estable. 



Por ejemplo, en el gráfico siguiente, tenemos dos puntos de equilibrio, x\ y 
X2, porque en ambos la fuerza es cero. Pero sólo el segundo es estable, porque si 
desplazamos una partícula desde X2 en cualquier dirección, la fuerza será hacia 
X2- Y por lo tanto habrá algún tipo de oscilación en torno a dicho equilibrio. 
Pero si una partícula la tenemos en x\, en reposo, y la desplazamos ligeramente, 
la fuerza será para alejar a la partícula de X\. 
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La diferencia entre ambos casos es que en torno a Xi la derivada ^f- < 

dF 
dx \x± 



dx \x 2 

(pendiente positiva), mientras que en torno a xi sucede lo contrario: 4ürl_ > 



(pendiente negativa). 

Veamos ahora qué sucede si en torno a un punto de equilibrio realizamos un 
pequeño desplazamiento Ax. 

En torno a x\ ya sabemos que no habrá oscilación. En torno a Xi si que 
la habrá, pero como la fuerza no es lineal con el desplazamiento, no será una 
oscilación armónica. ¿Cuánto vale la fuerza en el punto x = X2 + Ax? 

Desarrollo en serie de la fuerza F(x) en torno a un punto x 2 . 



F(x 2 + Ax) 



, N 1 dF 
^ + V.lx 



1 d 2 F 



Ax 2 + I.. 



entonces, si el desplazamiento Aa; es muy pequeño Aa; <C 1, se cumple que nos 
podemos quitar los términos Aa 2 , Ax 3 , etc, porque son mucho más pequeños 
que Aa;: Ax <C 1 => Ax 3> Ax 2 3> Ax 3 ^> . . . . Por lo tanto, en primer 
orden de aproximación, tenemos que: 



F(x 2 + Ax) 



dF 

dx 



Ax = — kAx 



como la derivada evaluada en x-i es una constante negativa, tenemos que la 
fuerza es lineal con el desplazamiento Ax. 

En resumen, cuando la fuerza F no es lineal podemos tener oscilaciones 
aproximadamente armónicas si: 

• se producen en torno a una posición de equilibrio estable 

• son el resultado de una pequeña perturbación con respecto al equilibrio 
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además, la constante recuperadora es la derivada de la fuerza evaluada en dicho 
punto de equilibrio estable. 

Podemos analizar también este tipo de oscilaciones en términos de energía 
potencial E p , en lugar de fuerza. Recordando la definición de energía potencial 
en fuerzas conservativas, tenemos que, si el sistema es unidimensional: 



F = — 



dE p 
df 



Siguiendo el caso de la fuerza representada en el gráfico anterior, tendríamos 
una energía potencial del tipo al que se muestra en la figura: 




Obsérvese que los máximos y mínimos relativos de la función E p (x) coinciden 
con los puntos donde la fuerza es cero. 



si x 
si x 



Xi 
%2 



dx 



.... . J inestable (E p (xi) max) 
[ estable (E p (^2) min) 



La diferencia es que mientras que en uno de los puntos de equilibrio, la energía 
potencial es un máximo relativo, en el otro ésta es un mínimo relativo. Si 
tenemos un máximo relativo (o local) de energía potencial, tenemos un equilibrio 
inestable, y si estamos en un mínimo relativo (o local) de la función E p {x) 
entonces tenemos un punto de equilibrio estable, en torno al cual puede haber 
oscilaciones. 

Desde el punto de vista de la energía potencial, para que dichas oscilaciones 
sean armónicas, la función E p (x) ha de ser una parábola del tipo 



E p = (l/2)fc(¡F - x cq f 



(3.7) 



En general esto no se produce, y en el caso de la gráfica tampoco. 
Vcramos qué sucede con la energía potencial si realizamos un pequeño de- 
splazamiento Ax desde la posición de equilibrio estable x<¿- Observar en la figura 
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que la referencia para la energía potencial se ha escogido como x rc f — x 2 , de 
modo que E p {x 2 ) = E p (x lc f) = 0. 

El desarrollo en serir de E p (x) en torno al punto X2 es: 



E p (x 2 +Ax) = E p (x 2 ) - 



1 dE p 
'Tí ~dx 



1 d 2 E p 

Ax> 



2! dx 2 



Ax 2 - 



1 d 3 E. n 



3! dx 3 



Ax 3 - 



1 
4! 



=0 



Obsérvese que los dos primeros términos desaparecen (en x 2 la energía po- 
tencial es cero, por ser la referencia escogida, y la primera derivada también 
lo es, puesto que se trata de un mínimo relativo). El primer término no nulo 
del desarrollo es el que va con Ax 2 . Si además suponemos el desplazamiento 
muy pequeño (Aa; <C 1) entonces nos podemos quedar con el primer término 
no nulo del desarrollo, porque Ax 2 3> Aa; 3 3> Ax A .... Entonces, podemos dar 
por buena la aproximación: 



A x 1 í d2E P 
E p{x + Ax),-(^-f 



Ax 2 



Con lo cual, identificamos la derivada segunda evaluada en el punto de equi- 
librio como la constante recuperadora k en la ecuación 3.7. Por tanto, la fre- 



cuencia es ui 



k/m 



(d 2 E p /dx 2 ) Xcq 



id 



3.8 Superposición de armónicos 

¿Qué sucede cuando dos oscilaciones armónicas se superponen? ¿Cuál es el 
resultado? Depende de cómo sean las superpuestas: 



3.8.1 Superposición de armónicos de igual frecuencia 

Tenemos dos armónicos de amplitudes y fases distintas, pero igual u>. El resul- 
tado de la superposición es otra oscilación armónica, con esa misma frecuencia: 



x\(t) = A\ sin(o;í + ai) 
x 2 (t) = A 2 sin(a;í + a 2 ) 



x(t) = Xi(t) + x 2 (t) = Asui(üjt + a) 



Veamos cuál es la amplitud resultante A, y la fase a. 

Utilizamos la relación trigonométrica de seno de la suma: sin(# + 
sin 6 eos ó + eos 8 sin ó. 



(3.8) 



x = A sin(cjí + a) = A (sin ujt eos a + sin a eos ojt) 
por otro lado: 

x = x\ + x 2 = Ai sin(cjí + a±) + A 2 sin(cjí + o¿ 2 ) => 
x = Ai (sin ojt eos ai + sinai coswí) + A 2 (sinwí cosa2 + sina2 coswí) (3.9) 
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Igualamos las ecuaciones (3.8) y (3.9), y sacamos por un lado factor común 
sinwí y por otro factor común coswí, es decir: 

sin Lot ( A eos a) + eos u>t(A sin a) = 

sin urt (A\ eos ai + Ai eos «2) + eos u>t (Ai sin a\ + Ai sin a-i) 

A continuación igualamos, por separado, la parte proporcional al seno y la 
parte proporcional al coseno. 



eos luí =>■ Asina = Ai sinai + A2 sina2 
sinwí =>■ Acosa = Ai cosai + A2 cosa2 
dividiendo la ecuación (3.10) entre la ecuación (3.11): 



(3.10) 
(3.11) 



tana 



Ai sin ai + A2 sin a2 
Ai eos ai + A 2 sin a2 



elevando al cuadrado las ecuaciones (3.10) y (3.11) y después sumándolas, se 
obtiene: 

A (eos a + sin a) = A 1 sin ai + A 2 sin a 2 + 2A1A2 sinai sina2+ 

+A\ eos 2 ai + A\ eos 2 a 2 + 2A X A 2 eos ai eos a 2 ==>■ 
A 2 = A\ + A\ + 2AiA2(cosai cosa2 + sin ai sina2) =£- 



A 2 = A 2 + A 2 2 + 2aiA 2 cos(a 2 - ai) 



Es interesante observar que la amplitud dependerá de la diferencia de fase 
entre las dos oscilaciones superpuestas. Si están en fase (ai = 012), entonces 
la amplitud es la suma A = Ai + A2, mientras que si están en contrafase 
(o¿2 — o¿\ = 7r), entonces la amplitud es la resta A = Ai — A2. 

3.8.2 Superposición de armónicos de igual amplitud y fase, 
pero distinta u 

Ahora vamos a estudiar un caso disinto: tenemos dos ondas de igual amplitud 
pero frecuencias distintas. La fase inicial la ponemos a cero por simplicidad. 



x\ = Asinwií 



x = x\ + X2 = A(sinu;ií + sinu^í) 



x 2 = Asin^íJ 
ahora echamos mano de la identidad trigonomética: 

sin 9 + sin ó — 2 sin ( — : — I o >-, 



esto nos permite escribir la ecuación previa así: 



x = A(sinu;ií + sinu>2t) = 2Acos 



t sin 
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Esta ecuación la podemos interpretar como una oscilación a frecuencia (wi - 
bj<¿)/2, (que es la frecuencia promedio) pero con una amplitud modulada, A ■■ 

2Acos[{ijj 1 -u 2 )t/'A- 




La amplitud thene unos valores que van desde hasta 2A. El periodo de la 
modulación de la amplitud es 2tt/(uji — lü 2 ). 

Se puede demostrar que si las amplitudes de las dos oscilaciones superpuestas 
son distintas, Ai ^ A 2 , entonces la amplitud toma valores en el intervalo Ai — 
A 2 < A < Ai + A 2 . 

Cuando las dos oscilaciones superpuestas están en fase, entonces A = Ai + 
A 2 — > Lüit = uj 2 t + 2n7r, con n número natural. 

Cuando las oscilaciones están en contrafase: A = Ai — A 2 , es decir: cuit = 
uj 2 t + (2n - l)ir. 



Chapter 4 

Sólido Rígido 



4.1 Introducción 

En este capítulo nos ocuparemos de la dinámica del sólido rígido, entendiendo 
como tal un sólido indeformable. Es decir, la distancia relativa entre las partículas 
que lo componen se puede considerar constante. 

El movimiento en general es una combinación de los dos tipos de movimiento 
simple: traslación y rotación. 

Para la traslación, podemos tratarlo teniendo en cuenta que, como vimos en 
los sistemas de partículas, su dinámica sería la de una partícula ubicada en el 
centro de masas. 

Sin embargo, si el sólido rígido rota en torno a un eje fijo, debemos introducir 
un nuevo elemento, el momento de inercia /. 



4.2 Rotación en torno a un eje fijo y momento 
de inercia 



4.2.1 Convenio de sentido de rotación 

Cuando un sólido rígido (o una partícula) rota con velocidad angular u) en torno 
a un eje, por ejemplo el OZ, existe un convenio según el cual se representa el 
vector lo — cok, paralelo al eje de rotación. El sentido de rotación se deduce 
utilizando la regla de la mano derecha, según la figura de a continuación: 

ÜJ 





55 
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donde el símbolo representa una flecha saliendo del plano donde está 
dibujado, apuntando hacia el lector, mientras que el símbolo (g) tiene la misma 
dirección pero sentido opuesto. 

Si escogemos un origen para el vector r¿ en el mismo eje, podemos escribir 
la velocidad de la partícula i como el producto vectorial: 

v l =uJxr l (4.1) 

Pensemos que |r¿ | sin 9 es el radio de la circunferencia que describe la partícula 
i, mientras que la velocidad w¿ es al mismo tiempo perpendicular a f¿ y a w. 
Todo tiene sentido pues. 



4.2.2 Momento angular de un sólido rígido en rotación 

Supongamos ahora que queremos averiguar el momento angular de una partícula 
de un sólido rígido que está rotando con velocidad angular ¡Jj || OZ, como se 
muestra en la figura. El origen de coordenadas O está ubicado en el eje de 

rotación. 

z 




Entonces, 

Jo.i = ñ><Pi = rmiñ x v^ (4.2) 

mientras que el momento angular total del sólido (que está compuesto por 
n partículas) sería: 

n n 

Jo = ^2 Jo,i = X! m i(n x Vi) 

[Recordar: Jq¿ es el momento angular de la partícula i (que forma parte 
del sólido rígido) con respecto al origen de coordenadas O. En cambio Jq es 
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el momento angular de todo el sólido rígido, es decir, la suma de los momentos 
angulares de cada una de las partículas, desde i = 1 hasta i = n] 
La proyección de Jo.i sobre el eje de rotación, que es el OZ, es: 

yo A = Jo,i ■ k 
Análogamente sucede con el momento angular total: 

í Jo) = Jo -k 
Esta última expresión la podemos desarrollar: 

n n n 

yJo) = ^2 i/ '*) = z2 m¿ ^ ' ^ x Vi) = ^2 mi ^ i ' $ x **») 

i— 1 í— 1 i— 1 

donde se ha utilizado para dar el último paso una propiedad de los productos 
escalar y vectorial: A- (B x C) = C ■ (A x B). 

Llegados a este punto podemos utilizar la expresión (4.1) para sustituir u¿ 
en la última ecuación: 

n n n 

Jo) = y~]miVi-(kxri) = ^m,((jxf,)'(/:xf 1 ) = ^m í w(fcxr¿)-(fcxr ¿ ) = 

i— 1 i—1 i—1 

n n / n \ 

= ^m ¿ w|£x fA 2 = ^m¿w(r¿ sin ¿ ) 2 = ^m ¿ <¿ 2 2 w 
¿=i ¿=i \¿=i / 

donde <¿¿ z = r¿ sin0¿ es la distancia de la partícula ¿ al eje de rotación. 

4.2.3 Momento de Inercia 

En resumen, vemos que tras unas pocas operaciones algebraicas podemos de- 
ducir que la proyección, sobre el eje de rotación, del momento angular es pro- 
porcional a la frecuencia lo y a una cantidad (entre paréntesis en la última 
expresión) que sólo depende de la estructura interna del sólido rígido. A esta 
cantidad la denominamos momento de inercia y la denotamos con la letra /. 





n 




U-- 


= E 

¿=i 


d 2 z mi 




w. ■ 


- I z v 



Observar que el momento de inercia depende del eje en torno al que gire el 
sólido. 

Si el sólido es continuo, es decir, es inviable realizar un sumatorio sobre todas 
sus partículas, entonces el momento de inercia se ha de calcular realizando una 
integral sobre todo el volumen del sólido. 



si n — > oo =>■ y mid i — > \ dmD 

¿=1 JVol 



2 
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donde la integral se realiza sobre todos los elementos diferenciales de masa dm 
que conforman el sólido, siendo D 2 el cuadrado de la distancia de cada uno de 
estos elementos al eje de rotación. 

En general no se puede escribir Jo — 13, porque el momento angular y el 
eje apuntan en direcciones distintas. Existen algunos ejes de rotación, que son 
ejes de simetría del sólido que está girando, tales que se cumple que J = IzCó. 
Dichos ejes también se suelen llamar ejes principales de inercia. 

Si queremos calcular el momento de inercia de un objeto bidimensional (por 
ejemplo un disco) la integral pasa de ser de volumen a ser una integral de 
superficie, mientras que si el objeto es lineal (por ejemplo una varilla) entonces 
la integral es una integral de línea. 



4.2.4 Momento de inercia de una varilla 



x = -L/2 



L/2 



Una varilla de sección despreciable, homogénea, de longitud L y masa M 
rotando en torno a su centro de masas: 

Antes de nada, la definición de momento de inercia, 



/ = / dmD 2 



cada elemento dm, ¿qué masa tiene? Los elementos infinitesimales de masa 
dm tienen todos la misma longitud, dx. 

Mediante una sencilla regla de tres (puesto que la varilla es homogénea) 
vemos que dm — Mdx/L. Vemos que estamos pasando de hacer una integral 
en la variable masa m a hacerla en la variable espacial x. Colocamos el origen 
de coordenadas en la intersección del eje con la varilla. La distancia de cada 
punto dm al eje es precisamente x. Finalmente quedan por definir los límites 
de integración: en nuestro caso la x va desde —L/2 hasta +L/2. Con lo cual: 



/ = / dmD 2 = 



L/2 , M 2 M 1 r 31 L/2 

dx — x = \x ó \ ' 



ML 2 
12 
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4.2.5 Momento de inercia de un disco 





Un disco plano, homogéneo, de masa M y radio R, rotando en torno a un 
eje que pasa por su centro de masas, y que es perpendicular al propio disco. 

Aquí tenemos un objeto bidimensional, con lo que la integral debe hacerse 
sobre una superficie. Lo más sencillo es utilizar coordenadas polares, y pasar 
de (x,y) a (r, 6), coordenadas que ya conocemos. La superficie total es wR 2 . 
Un elemento diferencial de superficie es, en polares (ver dibujo) rdrd6 1 y tendrá 
una masa: dm = rdrdOM /ttR 2 . La distancia D 2 es precisamente r 2 , y ahora las 
varaibles de integración son dos, r, 9. Sus límites de integración son < r < R 
y < 8 < 1-k. Entonces la integral del momento de inercia queda: 

I z = i dmD 2 = ¡^ de i" (r^r 2 ) dr = -^L [** ¿6 í" 
Jsu P Jo Jo V ^ 2 ) *R 2 Jo Jo 



r 3 dr 



M 
^R 2 



dO 



R i 



M i? 4 

^2tt — 

ttR 2 4 



MR 2 



4.2.6 Momento de inercia de un paralelepípedo 



h 
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Un bloque paralelepípedo, homogéneo, de masa M y dimensiones a—b—c, 
rotando en torno al eje que pasa por su centro de masas, y que está orientado 
en la dirección OZ. 

En esta ocasión tenemos un objeto tridimensional, la integral ha de hacerse 
sobre un volumen y por tanto en tres variables. En este caso las cordenadas 
cartesianas (x, y, z) son las más útiles. El volumen es V = abe, y la densidad 
p = M/V = M /abe. Un elemento diferencial de volumen sería dV — dx ■ dy ■ dz. 
Y la masa del mismo es: dm — pdV — Mdxdydz/abc. La distancia de cada 
punto (x,y,z) al eje de rotación es D 2 = x 2 + y 2 . Finalmente, los límites de 
integración son, si situamos el origen de cordenadas en el centro de la base: 
—a/2 < x < a/2, —6/2 < y < 6/2 y < z < c. La integral del momento de 
inercia nos queda así: 

r r c r h l 2 r a l 2 m 

I z = / dmD 2 = dz dy dx—-(x 2 +y 2 ) = 

Jvoi Jo J-b/2 J-a/2 abe 

M f c , f b ' 2 , fx 3 2 Y =a/2 M f c , f b l 2 , /a 3 2 

= ~T / dz / dy\—+xy\ = — / dz / dy — + ay 

abcj J_ b/2 V 3 J x =-a/2 abc h J-b/2 V 12 

M f c , fa 3 y ay 3 \ v=b/2 M f c , /a 3 6 a6 3 

dz —— H — — = — — / dz 



abe J \ 12 3 J _ b , 2 abe J \ 12 12 

M cab , 2 , 2 n M i 2 , 2 n 

= -; (a 2 + b 2 ) = — (a 2 + b 2 ) 

abe 12 y ' 12 V ; 

4.2.7 Teorema de Steiner 

Sea /* el momento de inercia de un sólido de masa M cuando rota en torno 
a un eje z que pasa por su centro de masas. ¿Qué sucede con el momento de 
inercia (1°) si rota en torno a un eje paralelo al mismo pero que no pasa por el 
centro de masas sino a una distancia D z del mismo? 

Según el teorema de Steiner, el momento de inercia cuando se gira en torno 
al nuevo eje es: 



i° = i: + md 



4.2.8 Ecuación del movimiento 

Si tenemos un sistema de n partículas, el momento angular Jo varía según la 
ecuación: 

d -§=Ñ = ±N f 

i 

como vimos en el tema II. Por otra parte la componente OZ de la ecuación 
anterior es: 

d Jo , , y . , 

notación 




%-»■ '«» 
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Por otra parte, en este tema hemos visto que para un sólido rígido rotando 
con frecuencia u = cok en torno a un eje fijo, la proyección del momento angular 
sobre el eje es: 



(Jo) 



= J z = I z u> 
z 



dJ z dtü 

-^=I z -r=I z a (4.4) 

dt dt K ' 

Combinando las ecuaciones (4.3) y (4.4) obtenemos una ecuación del movimiento: 



d 2 e 

Esta ecuación es el equivalente a la segunda ley de Newton para rotación de 
sólidos rígidos. A partir de la misma podemos deducir u>(t) y 9(t). 

Cuando el eje de rotación es uno de los ejes inerciales principales, y pasa por 
el centro de masas, se escribe con un superíndice * («/*, I*). Reescribimos la 
ecuación con que hemos dado comienzo a esta sección: 

di* 

dt 
Como estamos rotando en torno a un eje principal, se cumple que 

J* = ros 

Entonces, combinando estas dos ecuaciones y siguiendo un proceso análogo 
al anterior, deducimos la siguiente ecuación del movimiento: 



N* = la 



Esta ecuación es válida sólo cuando se rota en torno a un eje principal ele 
inercia. 

4.3 Energía Cinética de rotación 

La energía cinética de un sólido compuesto de n partículas, cada una de ellas 
con masa m¡ y velocidad w¿, es: 

1 n i n i n 

i i i 

puesto que consideramos que el sólido está rotando con frecuencia CS. Recorde- 
mos ahora la propiedad del producto escalar de un vector por el producto vecto- 
rial de otros dos: A-(B x C) = C ■ (Ax B) = B ■ (C x A). Entonces, retomando 
la última ecuación: 

E c = - ^2 mi ^ ' (^ x ^«) 

i 

como üj es igual para todos los puntos i, puede salir del sumatorio: 
Ec = o I 51 m% ^ i x Ví)>-lü= - J ■ CS 
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Por tanto en un sólido rígido rotando con ¿¡5, la energía cinética de rotación es: 

E c = ^J-io 



mientras que si tenemos rotación plana, J = Iü, y por tanto, sustituyendo 
esto en la ecuación anterior: 



1 9 

En = -Iu? 



4.4 Fuerzas de Ligadura 



Existen diversas situaciones en las que para mantener en el tiempo una rotación 
de un sólido rígido en torno a un eje fijo es necesario aplicar "fuerzas de liga- 
dura" . 



4.4.1 Caso I: Centro de masas no contenido en el eje de 
rotación 

Si el eje de rotación no pasa por el centro de masas i? c dm, entonces éste está 
rotando, y por lo tanto su aceleración no es cero. 

Ejemplo: Sólido rígido constituido por dos partículas de masa m situadas a 
distancia 2L. El eje de rotación es perpendicular a la línea que las une, pero no 
pasa por el centro de masas, sino que queda a una distancia L/2 del mismo. 

Lü 



III 





Fi 


edm 




L/2 





H F 2 



2L 

El centro de masas, por tanto, está describiendo un movimiento circular de 
radio L/2. El movimiento del edm es el equivalente al de una partícula de masa 
M = 2m que describe esa circunferencia, y su aceleración (ver final del tema 1) 
es A = oj 2 L/2. Por tanto, para mantener la rotación es necesario aplicar en el 
eje dos fuerzas de ligadura (F\ y F%) tales que F\ + F 2 = 2mui 2 L/2 = muí L, 
pero que al mismo tiempo no aporten momento. 

En definitiva, cuando el eje de rotación no pasa por el centro de masas, es 
necesario que actúen fuerzas externas para compensar la aceleración del centro 
de masas y que el eje de rotación se mantenga fijo. 
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4.4.2 Caso II: El eje de rotación no es un eje principal de 
inercia 

Si el eje de rotación pasa por el centro de masas pero no es un eje principal de 
inercia entonces suceden dos cosas: En primer lugar la aceleración del centro 
de masas es nula, puesto que el edm está fijo, y por el otro lado no se conserva 
el momento angular. Para preservar el giro, han de actuar una serie de fuerzas 
externas que hagan que se cumpla la ecuación de conservación del momento 
angular. 

Veamos un ejemplo utilizando el mismo sólido que en el caso I. 
I 




El centro de masas (O) no sufre movimiento de traslación alguno, con lo que 
su aceleración es cero, y por tanto no es necesario actúen fuerzas externas desde 
este punto de vista: j4 c( j m = — > "^2 F cxt = 0. 

Por otro lado el momento angular no se conserva (vemos que su dirección irá 
cambiando conforme vayan rotando las partículas de masa m que constituyen 
el sólido rígido). De hecho, para el observador situado en O, tenemos que Jo 
describe un movimiento circular, y por tanto: 



dJp 
~df 



lo x J 7^ 



Para que se mantenga la rotación del sólido rígido es necesario que actúe una 
serie de fuerzas externas tales que la mantengan. Es decir, se han de intro- 
ducir dos fuerzas externas de ligadura, i<\ y F 2 actuando sobre el eje, tales que 
el momento de las mismas cumpla la ecuación de conservación del momento 
angular: 



N» 



N F 



N Fo = 



dJo 
dt 



Pero al mismo tiempo, dichas fuerzas deben ser tales que F\ + F% = para 
que la suma de fuerzas externas no ocasione un movimiento de traslación del 
cdm. 
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4.5 Ejemplos 

4.5.1 El péndulo físico 

Anteriormente se estudió el péndulo ideal. Se trata de una idealización porque, 
entre otras cosas, se consdiera que la masa que oscila es puntual. En esta ocasión 
vamos a sustituirla por un péndulo más realista, con masa M pero también con 
un volumen finito que, además, puede ser irregular. El péndulo físico es un 
sistema que tiene un grado de libertad, que le permite rotar en torno a un eje. 




Supongamos que el sólido de la figura tiene momento de inercia I para el 
eje de rotación escogido. El punto fijo de rotación dista D del centro de masas 
(cdm). 

Por un lado tenemos la ecuación del movimiento: 

Ia = N 

¿Cuál es el momento N de la fuerza? En primer lugar, esa fuerza es el peso, 
como consecuencia del campo gravitatorio, y se aplica al centro de masas. El 
momento se define como el producto vectorial N = fx F. El vector fes el que va 
desde el origen escogido (en este caso el eje de rotación) hasta el punto donde 
la fuerza se aplica, es decir, es un vector de módulo D en nuestro diagrama. 
Entonces: 

N = -MgD sin9 

donde el signo negativo indica que la fuerza es recuperadora, en sentido 
opuesto a las variaciones angulares. 

Existiría otra fuerza, la fuerza de ligadura ejercida en el eje de rotación. Sin 
embargo el momento de esta fuerza es nulo, porque su distancia al origen es 0. 

Por otro lado tenemos la ecuación del movimiento: 



N = Ia = I 



dt 2 



De las dos últimas ecuaciones se deduce que: 

d 2 e 



dt 2 



-MgD sin t 
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Cuando el ángulo es pequeño, se puede aproximar por su seno: 

e <r i =*> sin e ~ e 

Por tanto, la ecuación del péndulo físico para pequeños ángulos es: 

d 2 9 MgD 

d¿- —° 

Es decir, para pequeños ángulos, el péndulo físico oscila armónicamente, con 
frecuencia w 2 = MgD/I, y un periodo T = lix^J 1 ¡MgD. 

4.5.2 Rodaduras sin deslizamiento 

La ecuación del movimiento de sólidos rígidos también nos permite resolver 
otros problemas complejos en los que intervengan sólidos rígidos rodando. Un 
ejemplo es el de una esfera que, rodando, desciende por un plano inclinado. La 
esfera tiene masa M y radio R. Como está rodando (y no desliza) está actuando 
en el punto de contacto con el plano una fuerza de rozamiento, dirigida contra 
el sentido de avance. Por otro lado, en ese mismo punto tenemos la reacción 
normal del plano, y además está el peso que se aplica en el centro de masas de 
la esfera. 

Con respecto a la fuerza de rozamiento, no es apropiado escribir una ex- 
presión del tipo F r — /j,N puesto que esa es la máxima fuerza de rozamiento. 
Sin embargo, como la esfera está rodando, no podemos afirmar que la fuerza sea 
la máxima (si estuviese deslizando sí que podríamos decirlo) 




Las fuerzas externas son, por tanto: 

F cxt = Mg, Ñ, F r 
Con los ejes que hemos escogido, descomponemos las fuerzas: 

Mg sin <pi 



Mg= . 

-Mg eos tpj 



N = Ni 

_£/ ip _£/ qp l. 



Tenemos por un lado un movimiento de traslación. La aceleración del centro 
de masas es A — Ai. Planteamos dos ecuaciones: 



E^: 



MA => i M9 /Z~ Fr = M n A ( 4 - 5 ) 

1 N — Mg eos (p = v ' 
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Por el otro lado, con respecto al eje de rotación (centro de masas), sabemos 
que sólo la fuerza de rozamiento aporta momento. 



N* 



N Fi =RF r (-k) 



El momento de inercia de una esfera maciza rotando sobre un eje que pasa 
por su centro de masas es: I* = (2/5)MR 2 , y la ecuación del movimiento de 
rotación de los sólidos rígidos nos dice que N* — I* a, con lo que, utilizando la 
ecuación anterior: 



N* = RF r (-k) 
Ñ* = 7*5 



Fr 



: MRa 



(4.6) 



puesto que a = a(—k). 

Tenemos de momento tres incógnitas, la fuerza de rozamiento, la aceleración 
angular a y la aceleración del centro de masas. Pero sólo tenemos dos ecua- 
ciones que nos sean útiles (porque la segunda ecuación de (4.5) no nos aporta 
información). Una última ecuación la planteamos teniendo en cuenta que la 
rodadura es sin deslizamiento, por tanto 



A = aR 



(4.7) 



Ahora ya tenemos tres ecuaciones (4.5), (4.6) y (4.7) y tres incógnitas. 

Resolviendo el sistema se obtiene que A = (5/7)gsimp, a = (5/7)gsmip/R 
y F r = (2/7)Mgsinip. 

Hemos supuesto que hay rodadura sin deslizamiento, es decir, que la fuerza 
de rozamiento es menor que la fuerza de rozamiento máxima: 



7 



Mg sin ip < ¡j,N — \xMg eos ip 



Es decir, habrá rodadura siempre y cuando el coeficiente de rozamiento sea 
¡i > = tan</?. 

Rodadura sin deslizamiento: relación de A y a 

Supongamos que tenemos una esfera, disco o cilindro rodando sin deslizar por 
un plano, tal y como se muestra en la figura. 
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Si el sólido rígido rota un ángulo d9, vemos que el punto P avanza una 
distancia que es la misma que el arco que ha rodado. El arco es dX = rdO, 
donde r es el radio (importante! el ángulo en radianes). Entonces: 

dX = rd6 

Derivando con respecto al tiempo sacamos la velocidad del punto P. 



dX 



de 



dt r dt 



V = ruj 



Si volvemos a derivar, tenemos la aceleración: 



d 2 X 
dt 



dt 2 



A = ra 



4.5.3 La máquina de Atwood 

La máquina de Atwood consta de dos bloques de masas mi y mi < mi , unidas 
por una cuerda a través de un disco de masa M y radio R que hace las veces de 
polea. 




m x g 



Suponemos que la cuerda, inextensible y sin masa apreciable, no desliza. 

Antes de plantear ecuaciones, notar que, debido a la diferencia de masas, las 
aceleraciones de los bloques 1 y 2 son iguales y de sentido contrario. 

En primer lugar, planteamos la 2 a ecuación de Newton para los bloques de 
masa mi y mi: 

m\g — Ti = mía (4-8) 



Tn 



mig = mía 



(4.9) 



Por otro lado, la polea está rotando en torno a su eje, con w \\ k. Sobre la polca 
actúan varias fuerzas: su propio peso, la fuerza de ligadura que la sostiene, y las 
tensiones Ti y Ti. Desde su centro de masas, sólo las dos últimas contribuyen 
a su momento N* . 

N* =TiR-T 2 R 

Por otro lado, el momento de inercia de un disco de masa M y radio R es 
I* = MR 2 /2, y la ecuación del movimiento de los sólidos en rotación es: 

N* = I*a = ]-MR 2 a 
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Juntando estas dos ecuaciones: 

(Ti - T 2 )R = -MR 2 a (4.10) 

De momento tenemos las incógnitas de las tensiones Ti y T2, la aceleración o 
de los bloques y la aceleración angular a de la polea. Necesitamos una ecuación 
más. 

Supongamos ahora que estamos en un punto del borde del disco. Como la 
cuerda rueda sin deslizamiento, planteamos la ecuación de rodadura: 

a = aR (4.11) 

La aceleración en el borde del disco es el producto del radio por la aceleración 
angular. 

Ahora ya tenemos 4 ecuaciones: las (4.8), (4.9), (4.10) y (4.11). Resolvemos: 

(4.8) + (4.9) -»• 

(mi — m 2 )g + {m\ + m 2 )a = T\ — T 2 
Con las ecuaciones (4.10) y (4.11) obtenemos la expresión: 

(Ti - T 2 ) = l -Ma =► 
sustituyendo el la ecuación anterior: 

(mi — rn,2)g + {mi + m 2 )a = -Ma =*> 

mi — m 2 
a = g- 



TOi + m 2 + M/2 

Es la aceleración con que la masa mayor cae y la masa menor sube. La polea 
rota con una aceleración angular: 



mi — m 2 
a = g- 



R{mi +m 2 + M/2) 



Chapter 5 

Termodinámica 

5.1 Introducción 

El término termodinámica procede de las voces griegas "termo" (calor) y "di- 
namo" (fuerza) y, literalmente hace referencia a los procesos mecánicos inducidos 
por el calor, aunque el campo de la termodinámica es bastante más amplio. Así, 
sirve para estudiar el funcionamiento y rendimiento mecánico de máquinas, pero 
también de sistemas químicos, biológicos, geológicos. 

La termodinámica es una ciencia experimental en la que se estudian sis- 
temas macroscópicos. Un sistema puede ser cualquier porción del universo, 
esté o no aislada del entorno. Para estudiar y caracterizar un sistema ter- 
modinámico se echa mano de una serie de magnitudes que se evalúan sobre sis- 
temas macroscópicos. De esta manera se dice que las variables termodinámicas 
caracterizan el estado de un sistema. Dicho de otro modo, las variables ter- 
modinámicas son funciones de estado y que están en equilibrio termodinámico. 
Las variables termodinámicas solo pueden definirse cuando el sistema se encuen- 
tra en equilibrio termodinámico. Esto comporta tres equilibrios simultáneos: 
equilibrio térmico, mecánico y químico. 

Los sistemas termodinámicos pueden experimentar cambios (por ejemplo de- 
bido a un aporte de calor, o a la realización de un trabajo...). A estos cambios 
los denominaremos procesos termodinámicos, y pueden ser de muy diversa natu- 
raleza, pero están todos regidos por los principios de la termodinámica, algunos 
de los cuales veremos a lo largo de este capítulo. 

5.2 Equilibrio térmico y temperatura 

El primer concepto que nos evoca el término termodinámica es el de tempera- 
tura. La temperatura es la magnitud física que relacionamos con los conceptos 
comunes de "caliente" o "frío". 

Decimos que dos sistemas están en equilibrio térmico si se encuentran a la 
misma temperatura. 

Principio Cero de la Termodinámica Dos sistemas en equilibrio térmico 
con un tercero están en equilibrio térmico entre sí. 

69 
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Así se puede definir una escala de termperaturas, utilizando un termómetro, 
es decir, un sistema que posea alguna propiedad que dependa de la temperatura, 
y que a la vez pueda ser mcdible. El ejemplo más común de termómetro es el 
de mercurio. Su dilatación con la temperatura nos permite definir una escala: 
La longitud de la columna de mercurio varía con la temperatura. 

Existen varias escalas, que pueden ser definidas tomando dos puntos fijos a 
los que se asigna una temperatura: 

• Escala Celsius: Se toma el punto (o temperatura) de congelación del agua 
como o C, y el punto de ebullición del agua como 100° C. Si disponemos 
de una columna de mercurio, marcamos la altura que ésta tiene a los o y 
100° C y el espacio entre las dos marcas se divide en 100 intervalos iguales. 
De esta manera se define la escala. Es la más empleada en situaciones 
cotidianas. 

• Escala Fahrenhcit: Los puntos de referencia son distintos. El de fusión 
del hielo es 32° F y el de ebullición del agua es 212° F. Si la temperatura 
en escala Fahrenheit la escribimos como Tp y la temperatura en escala 
Celsius como Te, ambas se relacionan así: 

T c = l(T F -32) 

Es empleada fundamentalmente en países anglosajones. 

• Escala Kelvin: También denominada escala absoluta. Es la empleada en 
el sistema intenacional. Lord Kelvin estableció que la temperatura límite 
o cero absoluto es —273.15° C. A partir de ahí se define la temperatura 
en Kelvin T¡<: 

T K = -273.15 + T c 



5.3 Leyes de los gases ideales 

Las leyes de los denominados gases ideales responden a una serie de observa- 
ciones empíricas históricas: 

5.3.1 Ley de Boyle/Mariotte 

A temperatura constante, y para gases a baja densidad, se encuentra que el 
producto entre la presión P y el volumen V es: 

PV = cíe 

Este comportamiento fue observado por Boyle (1662) y por Mariotte (1676) 
en estudios independientes. 



5.3.2 Ley de Gay-Lussac/Charles 



A presión constante, y, siempre en condiciones de baja densidad, se encuentra 

que el cociente entre el volumen ocupado por un gas V y su temperatura T 

también es constante: 

V 

j, = cíe 
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Se suele atribuir este descubrimiento a Gay-Lussac (1802) aunque ya lo había 
estudiado Charles (1787). 

5.3.3 Avogadro y el mol 

Una vez aceptada la hipótesis atómica, se introduce el concepto de mol. Un mol 
de cualquier sustancia contiene el mismo número de átomos (o moléculas). Este 
número es el conocido como número de Avogadro: 

N A = 6.022 • 10 23 partículas/mol 

Un mol de cualquier sustancia contiene siempre iV¿ átomos (o moléculas) 
y su masa viene caracterizada por el número másico del núcleo del elemento o 
compuesto en cuestión. 

5.3.4 Ley de Avogadro 

Avogadro encontró que a presión y temperatura constantes, el cociente entre el 
volumen ocupado V y el número de moles n es también constante: 

V 

— = cíe 
n 

5.3.5 Ecuación de Estado 

Las anteriores leyes se resumen en la siguente ecuación de estado propuesta por 
Clapeyron en 1834: 

PV = nRT 

donde P es la presión, V el volumen, n el número de moles y T la tem- 
peratura. La constante de proporcionalidad R es la constante universal de los 
gases 

i? = 8, 314 J/molK 

R = 0.082 atm-1/mol-K 

5.3.6 Ley de Dalton 

Es habitual encontrarse en situaciones donde en un recipiente existe no un gas 
puro sino una mezcla de varios conviviendo simulatáneamentc. Por ejemplo, 
si se encierra en un cilindro aire procedente de la atmósfera (a la altura del 
nivel del mar) tendríamos una mezcla de nitrógeno y oxígeno fundamentalmente, 
además de pequeñas proporciones de otros gases como el dióxido de carbono, etc. 
Entonces, suponiendo que se dan las condiciones de gas ideal, ¿cómo utilizar la 
ecuación de los gases ideales? En primer lugar definimos el concepto de presión 
parcial como aquélla que ejercería cada uno de los gases de la mezcla si éste 
estuviese solo ocupando el volumen total y a la misma temperatura. 

Si tenemos un recipiente de volumen V a temperatura T con NG gases, la 
presión parcial del 'i'-ésimo de ellos es P¿. Si hay n¿ moles del gas 'i'-ésimo, 
entonces 

PV = n,RT 



Fundamentos de Física I - Guillermo Sánchez 

La presión total Ptot es la suma de las presiones parciales, 

NG NG 

V 



E^-^v RT 



La ley que establece que la presión de una mezcla de gases es la suma de las 
presiones parciales de los componentes de la mezcla se denomina ley de Dalton. 
Es válida para mezclas de gases a la misma temperatura y que no reaccionan 
químicamente entre sí. 



5.4 Teoría cinética 

Consideremos ahora a un gas, que ocupa un volumen V, que consta de N 
moléculas. Vamos a realizar varias hipótesis: 

1. Las moléculas son puntuales, y tienen una velocidad v — (v x ,v y ,v z ). 

2. No existe interacción a distancia entre las partículas. 

3. Colisionan entre sí elásticamente. 

4. No existen direcciones privilegiadas: (v 2 .) = (v 2 ) = (v 2 ) = (1/3) (v 2 ). 



5.4.1 Energía cinética 

En función de la temperatura la velocidad de las partículas tendrá una dis- 
tribución u otra. Boltzmann estudió estas distribuciones para demostrar que 
estadísticamente la velocidad cuadrática media (v 2 ) se relaciona con la tempe- 
ratura: 

2 3k B T 

(v ) = 

m mo i 

donde T es la temperatura, m mo i la masa de las moléculas y k B es la cons- 
tante de Boltzmann y vale k B = 1,381 • 10~ 23 J/K. 

Entonces, si tenemos n moles de gas, y por lo tanto uNa partículas, la 
energía cinética es: 

1 1 Q / T Q 

E c = nN A -m mo \(v 2 ) = nN A-ro^sX = nN A -k B T 

2 2 jrinisr 2 

Si escribimos Na ■ k B — R tenemos que 

E c = n ■ 3 RT 

2 

En su movimiento, las partículas tienen 3 grados de libertad. Cada grado 
de libertad aporta la misma proporción de energía cinética (porque no hay 
direcciones privilegiadas). Por tanto, cada grado de libertad contribuye en un 
factor (l/2)k B T. 
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5.4.2 Presión 

Por otro lado, veamos cómo es la presión en los gases ideales. La presión se 
define como fuerza dividida por superficie. 

Tomemos un pequeño lapso de tiempo Ai. A partir de aquí podemos definir 
un elemento de volumen como el observado en la figura: un paralelepípedo de 
sección A y con las aristas perpendiculares a dicha sección miden v x At, que es 
la distancia que una partícula del gas con velocidad v x puede recorrer en ese 
intervalo de tiempo. 




-A 



¿Cuál es la presión que el gas ejerce sobre una de las caras de ese cubo? 
Tomemos por ejemplo el área marcada (A) en la figura. La presión será la 
componente de la fuerza en la dirección perpendicular a esa superficie (dirección 
x) dividida precisamente por la superficie: 



A 



1 dp x 
A dt 



1 Ap x 
A Ai 



donde hemos tenido en cuenta que la fuerza es la derivada con el tiempo del 
momento lineal, para después aproximar la derivada por la variación con el 
tiempo, suponiendo que Ai es pequeño. ¿Cómo varía el momento lineal? Si 
suponemos que cada partícula, de masa m, tiene una velocidad v x y choca 
con otra de misma velocidad y sentido contrario, tras la colisión elástica su 
velocidad pasará a ser — v x . Es decir, su velocidad ha variado en 2v x . Y por 
tanto su momento ha variado Áp x = 2v x m. 

Si denominamos N' al número de partículas que en el intervalo de tiempo 
Ai chocan en A, entonces, el cambio de momento lineal de todas las partículas 
es 

Ap x = N'2v x m 



y por tanto la presión es: 



P 



1 Apx 
A Ai 



1 2mv 



A Ai 



-N' 



N es el número de partículas hay en el sistema de volumen V . El elemento 
de volumen que estamos estudiando tiene un volumen V = v x At ■ A. El número 
de partículas que hay en nuestro elemento de volumen se obtiene sencillamente 
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mediante una regla de tres. Pero solo la mitad de éstas partículas tendrán 
velocidad en la dirección x > 0, y por tanto sólo la mitad de esas partículas 
llegarán a la superficie A (las demás se van en sentido contrario, alejándose). 
Ese es el número N' que buscamos: 

N 1 

N' = —Av x At- 

V 2 

Entonces sustituimos este valor en la fórmula de la presión: 



1 %mv x N y 1 N 2 



Como estamos estudiando un número muy elevado de partículas, es apro- 
piado trabajar con cantidades promediadas: En lugar de v x trabajaremos con 
{vi). 

Según la hipótesis de isotropía (ausencia de direcciones privilegiadas) (v^) = 
(l/3)(w 2 ). Y, según hemos visto en el estudio de la energía cinética de los gases, 
el promedio de la velocidad cuadrática es 

2 3k B T 

{v = 

m 

con lo cual, la presión la podemos reescribir como sigue: 

N 1 tk B T N , 

Recordar que el número total de partículas es igual al número de Avogadro 
Na multiplicado por el numero de moles n, es decir, 

P = Kb-Í 

V 

que es la misma ecuación que la de los gases ideales si definimos R = 
N A k B => 



PV = nN A -k B -T = n • RT 

Es decir, mediante el análisis estadístico, la teoría cinética (que usa las 
hipótesis detalladas al inicio de este punto) nos permite desarrollar una ecuación 
equivalente a la de los gases ideales. 



5.4.3 Energía Interna 

En el modelo que obtenemos mediante el estudio cinético de los gases, se está 
suponiendo que la energía que tiene nuestro sistema es la energía cinética de las 
partículas que lo componen. La velocidad de las partículas depende solamente 
de la temperatura. Por tanto, la energía interna se puede escribir como función 
deT. 
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5.5 Energía Interna, Calor y Trabajo: Primer 
Principio de la Termodinámica 

Antes hemos visto que la energía de un gas ideal está asociada a su temperatura. 
Se trata de la energía interna y la denotamos con la letra U mayúscula. 

U = \nRT 

Si el sistema está aislado, su energía interna permanecerá constante. Pero si 
no lo está pueden producirse cambios en la misma. 

5.5.1 Calor 

Entre dos sistemas termodinámicos a distinta temperatura y conectados se pro- 
duce una transferencia de energía. En concreto hay una transferencia de energía 
procedente del sistema que se encuentra a mayor temperatura hacia el de tempe- 
ratura más baja. La transferencia existe hasta que ambos sistemas se encuentran 
en equilibrio termodinámico entre sí. Denominamos calor Q al flujo de energía 
como consecuencia de una diferencia de temperaturas. 

Al tratarse de una magnitud de energía se mide en Julios. Otra unidad 
empleada es la caloría: 

1 cal. = 4.184 J 

La caloría es la cantidad de calor que es necesario aportar a un gramo de 
agua para incrementar su temperatura en un grado Cclsius o Kelvin. 

5.5.2 Cambios de temperatura 

De manera similar a la definición de caloría, se puede definir la capacidad 
calorífica C de cualquier sistema como la cantidad de calor que es necesario 
aportar para incrementar su temperatura en un grado. Del mismo modo, la 
cantidad de calor para provocar un incremento de la temperatura del sistema 
en AT sería: 

Q = CAT 

en unidades del SI C se mide en J/K. En ocasiones se emplean candidades 
por unidad de masa o por mol. En ese caso se habla de capacidad calorífica 
específica y se utiliza una letra minúscula c. 

Cabe destacar que el calor Q es una magnitud escalar que lleva un signo aso- 
ciado a la ganancia o pérdida de energía interna. El sistema cuya temperatura 
disminuye está cediendo calor, y por lo tanto Q < 0, mientras que el sistema a 
menor temperatura está recibiendo calor, es decir, Q > 0. 

La expresión anterior para el calor está tomada considerando que C (o c) no 
cambia con la temperatura, es constante. En muchos materiales esto se cumple 
para amplios rangos de temperaturas, pero no es una ley general. Si C (o c) 
cambia con T, entonces lo adecuado es escribir el flujo de calor para un cambio 
infmictesimal de temperatura dT: es decir, el sistema pasa de una temperatura 
T a T + dT. El calor recibido es: 

dQ = CdT — racdT 

donde m es la masa. 
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5.5.3 Cambios de Fase 

Los cambios de fase son procesos termodinámicos en los que no hay cambios de 
temperatura pero sí un aporte (o pérdida) de energía en forma de calor. Un 
bloque de hielo a o C en una habitación a 20° C está recibiendo un aporte 
de calor que está haciendo incrementar su energía interna. Si el aporte ener- 
gía se mantiene se producirá una transición de fase: el hielo se funde. Otros 
ejemplos son la solidificación, la vaporización, la condensación... . Durante las 
transiciones de fase la temperatura del sistema no varía porque toda la energía 
recibida por el sistema es empleada en, por ejemplo, romper los enlaces del 
estado sólido del hielo. El calor que es necesario para provocar el cambio de 
fase de fusión de un bloque de masa m es: 

Qf = mL f 

donde L¡ es el denominado calor latente de fusión. El calor latente por mol se 
escribe lf, y entonces, para fundir n moles hace falta Qf — nl¡. Análogamente 
sucede con los calores latentes de vaporización, condensación... . 

5.5.4 Trabajo 

Los sistemas termodinámicos también pueden cambiar su energía interna rea- 
lizando (o absorbiendo) trabajo. Supongamos que tenemos un sistema ter- 
modinámico encerrado en un cilindro con un pistón de área A. 

Si el gas realiza una fuerza sobre el exterior que provoca un desplazamiento 
dx del pistón, entonces está realizando trabajo. 

dx 




Por la definición de trabajo 

W = F -dr= I F T dx = / PAdx = I PdV 



El trabajo es la integral de la presión en volumen. El elemento diferencial 
de trabajo es: 

dW = PdV 

Cuando el sistema realiza un trabajo hacia el exterior, el criterio que emplea- 
mos es el de W > 0. Si, en cambio, el sistema recibe trabajo desde el entorno, 
entonces decimos que éste es negativo W < 0. 

Nota: Existe otro criterio de signos, según el cual el trabajo efectuado sobre el 
sistema es positivo y negativo el trabajo que el sistema realiza sobre el exterior. 
En ese caso el elemento diferencial de trabajo es dW = —PdV. Sin embargo 
nosotros emplearemos el otro criterio de signos, enunciado anteriormente. 
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5.5.5 Primer Principio de la Termodinámica 

La variación de la energía interna de un sistema (AC7) es igual al calor Q trans- 
ferido mas el trabajo W realizado por (o sobre) el sistema, es decir: 



AU=Q- W 



En termodinámica se estudian los procesos llamados "cuasiestáticos" : proce- 
sos en los cuales el sistema evoluciona atravesando siempre estados de equilibrio. 
Así, todos los puntos del sistema tienen la misma presión, temperatura... . 

Así, por ejemplo a la hora de calcular el trabajo en un proceso de ex- 
pansión, el elemento diferencial de trabajo lo calculamos considerando que a 
cada pequeña variación de volumen dV, el sistema alcanza el equilibrio, y des- 
pués experimenta la siguente pequeña variación dV, etc.. hasta alcanzar el 
volumen final. 

En procesos cuasiestáticos, suficicntcmcnt lentos de forma que nunca se 
abandona el equilibrio termodinámico, resulta conveniente la formulación di- 
ferencial: 

dU = dQ- dW 

5.5.6 Experimento de Joule 

En 1844 Joule propuso un experimento esquematizado en la figura: Dentro 
de un recipiente que no permite el flujo de calor (un calorímetro) se ubican dos 
recipientes Ay B conectados mediante una válvula R, además de un termómetro 
para medir cambios de temperatura. Si A contiene un gas ideal y B está vacío, 
entonces al abrir la válvula R el gas ideal se expandirá para ocupar todo el 
volumen. 




Termómetro 



Figure 5.1: Esquema del experimento realizado por Joule en 1844 

El sistema está aislado adiabáticamente, es decir, no permite el flujo de calor. 
Por otra parte, el sistema no está realizando ningún trabajo hacia el exterior ni 
tampoco lo está recibiendo. Por tanto Q — W = 0. Entonces, según el primer 
principio ¿SU = 0. 

Joule encontró que en este proceso en el que cambian la presión y el volumen 
del gas ideal, su temperatura permanece constante. Por lo tanto dedujo que la 
energía interna de un gas ideal solo depende la temperatura, tal y como habíamos 
predicho anteriormente. 
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5.6 Procesos Termodinámicos 

Es costumbre representar tanto los estados termodinámicos como los procesos 
experimentados por el sistema según diagramas PV (Presión- Volumen). Estos 
esquemas facilitan el análisis de algunos procesos fundamentales que enumera- 
mos a continuación: 

1. Isóbaros: Presión constante. 

2. Isócoros: Volumen constante. 

3. Isotermos: Temperatura constante. 

4. Adiabáticos: Sin flujo de calor, Q = 0: 

A continuación comentamos brevemente algunas características de estos pro- 
cesos: 

5.6.1 Procesos Isócoros 

Los procesos isócoros o isócoros son aquellos en los que el volumen permanece 
constante. 



Po 



P 
Pi 



V 



Figure 5.2: Diagrama P-V: Ejemplo de proceso isocoro 

En los procesos isócoros, al no haber variación de volumen, el trabajo es nulo. 
Por tanto, toda variación de la energía interna del sistema será consecuencia del 
flujo de calor 

W = =*> AU = Q 

La capacidad calorífica de una sustancia, en general, se define como la 
variación de calor bajo variaciones de temperatura 



C 



dQ 

dT 



y puede depender tanto de la naturaleza del sistema como de las condiciones P, 
T,V). 

Cuando trabajamos a volumen constante se define C v , la capacidad calorífica 
a volumen constante. 

Para los gases ideales la energía interna solo depende de la temperatura. En 
estos casos se encuentra que el diferencial de energía interna es directamente 
proporcional al diferencial de temperatura: 



dU = dQ = C v dT 
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donde C v es la denominada capacidad calorífica a volumen constante. A veces 
es más cómodo utilizar la capacidad específica c v — C v /n. Depende de cada 
sustancia y, en principio, puede variar también con la temperatura y la presión. 
No obstante, en multitud de procesos es aproximadamente constante, con lo que 
se puede integrar de manera sencilla la expresión anterior: 

AU = Q = C V AT = nc v AT 

Por otra parte, si trabajamos con un gas ideal a temperatura inicial Ti, 
entonces la variación de temperatura cuando se pasa de una presión P\ a una 
presión P 2 es: 

PoV PiV V 

AT = T 2 -T 1 = -^— - — = <P 2 - Ti) — 
nR nR v ' nR 

5.6.2 Procesos Isóbaros 

Al ser procesos en los que la presión permanece constante, el trabajo total se 
calcula de forma inmediata. En un proceso isóbaro en el que el volumen inicial 
sea Vi y el final V 2 , a presión P, el trabajo es: 

W=í PdV = p[ dV = P(V 2 -V 1 ) = PAV 
Jv, Jv, 



P 



Vi V V 2 

Figure 5.3: Diagrama P-V: Ejemplo de proceso isóbaro 

Si el proceso isóbaro está protagonizado por un gas ideal, entonces también 
podemos tener en cuenta la ecuación de los gases ideales para calcular cual va 
a ser la variación de temperatura del sistema: 

PV 

PV = nRT => T = 

nR 

Si el sistema está cerrado (no entra ni sale materia y n es constante), y la 
temperatura inicial es Ti, entonces la variación de temperatura es: 

PVo PVi P 

AT = T 2 -T 1 = — - = {V 2 -V 1 ) — 

nR nR nR 

Cuando se trabaja a presión constante, se define la capacidad calorífica a 
presión constante o C p . También es habitual trabajar con la capacidad calorífica 
especifica a presión constante, c p minúscula que es la anterior dividida por el 
número de moles: C p = nc p . 
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Entonces los incrementos diferenciales de calor son: 

dQ = CpdT = nCpdT 

si consideramos que la capacidad calorífica a presión constante no varía con 
la temperatura, entonces: 

Q = C P AT = nc p AT 

Relación entre C v y C p para gases ideales 

En una transformación a volumen constante: 

dU = dQ = nc v dT 

En una transformación a presión constante: 

dU = dQ - dW = nc p dT - PdV 

La variación de energía interna no depende del tipo de transformación. La 
energía interna es propiedad del sistema (y no del proceso seguido), por lo tanto 
la variación de energía interna solo depende del estado inicial y final. Entonces, 
se pueden igualar las dos ecuaciones anteriores: 

nc v dT = ncpdT — PdV 

En un gas ideal, por la ecuación de estado, sabemos que a presión constante 

PV = nRT =¡> PdV = nRdT 

sustituimos en la expresión previa: 

nc v dT — nCpdT — nRdT =>■ nc v = nc p — nR 

Es decir, en los gases ideales se cumple que: 

nc v = ncp — nR => C' v = C p — nR 

Además, en gases monoatómicos y diatómicos se pueden calcular los valores 
teóricos de las constantes. 

• Gases monoatómicos: C v — |ni?, C p = |ni?. 

• Gases diatómicos: C v = |ni?, C p = |ni?. 

Por otra parte se define el coeficiente adiabático 7 = -J-. 

5.6.3 Procesos Isotermos 

En los procesos isotermos, la temperatura es constante. Que la temperatura sea 
constante no implica la ausencia de transferencia de calor. 

Si el gas es ideal, entonces PV = constante, por la ecuación de los gases 
ideales. 
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Vi v v 2 

Figure 5.4: Diagrama P-V: Ejemplo de proceso isotermo 



La energía interna de los gases ideales depende solo de su temperatura. Por 
tanto, al no haber variación térmica => AU = 0. Por tanto, 

Q = W 



Como bien es sabido, PV = k (de), es decir, 
PV = nRT 



p _ nRT 

r V 



w=OPdv 



w 



V2 nRT Í V2 dV 

dV^nRTl ^L^nRT[\nV} 



Vi 



V 



V 



V 



W = nRT(\nV 2 -\aVi) = nRT ln 



V2 
Vi 



5.6.4 Procesos Adiabáticos 

En los procesos adiabáticos no existe transferencia de calor. Por tanto toda 
variación de energía interna será consecuencia de la realización de un trabajo. 




Figure 5.5: Diagrama P-V: Ejemplo de proceso adiabático 

Análogamente al caso anterior: la no transferencia de calor no implica con- 
servación de la temperatura. 

Sabemos que dQ = y que Q = 0. 



dU 



-pdV 



Si el proceso adiabático lo experiment un gas ideal, como vimos anterior- 
mente, siempre podemos escribir la variación de energía interna en función de 
c v . Por otra parte, de la ecuación de los gases ideales, podemos despejar la 
presión P. 
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d Z = ~ Pd Z 1 , dV dT RdV 

dU = nc v dT \ => nc v dT = -nRT — =*> — = 

PV = nRT =► P=^J F T Cv V 

Esta expresión se puede integrar. El resultado nos da sendos logaritmos a 
cada lado del igual: 

lnT= InV + cte =$■ 

c v 

lnT + lnF^ =cte =*• 

TV^ = cte 

Además, por la ecuación de los gases, T = —5-, y sustituyendo en la ecuación 
anterior: 

PV^ +Í = cte 

Como se vio anteriormente, c p — c v + R — > 

_ + l = Cp ~ C " +1 = °p = 7 

C^j C^, C^> 

donde 7 es el coeficiente adiabático. Por tanto, podemos escribir que en los 
procesos adiabáticos con gases ideales se cumple que: 

PV 1 = cte (y, además) TV 1 ^ 1 = cte 

La primera de estas expresiones nos servirá para calcular el trabajo realizado 
por un gas ideal en procesos adiabáticos. Si el volumen que ocupa inicialmente 
el sistema es V\ y al final V^: 

w = Í V2 Pdv = i 2 c -^dv = -*t- (vr +1 vr +1 

Llegados a este punto, recordemos que, si la presión y temperatura iniciales 
son P\ y T\ respectivamente, y las finales son P 2 y T 2 , entonces: 



W 


= PiVl = cte 


PiV? ■ V{ 


^ +1 = cte ■ Vp +1 


PxVx 


= cte ■ Vp +1 



P-íY-í = nRTx =*> nRTx = cte ■ Vp +1 
Y análogamente con P 2 ^2- Entonces el trabajo queda: 

w = ^TT ( v '" +1 - v ^ 1 ) = =^Vt «K - **> = =7$T {T > - Tí > 

Finalmente, por la definición de 7 y la relación entre las capacidades caloríficas 
específicas, obtenemos: 

W = -nc v (T 2 - T^ 
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5.7 Máquinas Térmicas y Segundo Principio de 
la Termodinámica 

En termodinámica se pueden dar dos tipos de procesos si nos atenemos a su 
reversibilidad: Los procesos irreversibles y los reversibles. Los primeros son 
procesos en los que, una vez alcanzado el estado final, nunca se podrá regresar 
(de manera natural) al estado inicial. Por contra, en los reversibles, el estado 
inicial se puede alcanzar nuevamente. En los procesos irreversibles se observa 
que ha aumentado el desorden. La magnitud termodinámica que sirve de me- 
dida cuantitativa del desorden de un sistema es la entropía (S), de la que nos 
ocuparemos más adelante. 

5.7.1 Procesos espontáneos 

Hasta ahora no nos hemos ocupado de si un posible proceso termodinámico es 
espontáneo (es decir, se produce de manera natural) o no. Nuestra experien- 
cia nos sirve muchas veces para considerar si un proceso se puede dar o no de 
manera espontánea. Por ejemplo sabemos que una rueda girando, puede dejar 
de hacerlo por efecto del rozamiento tras un tiempo, aumentando así su tempe- 
ratura (energía cinética— > Calor), pero el proceso opuesto no sucederá: Nunca 
se enfriará para producir energía cinética, empezando a girar espontáneamente 
desde el reposo. El proceso anterior era irreversible. 

Análogamente sucede con la expansión libre de un gas (ver experimento de 
Joule). Si un gas está encerrado en una cápsula, comunicada con una válvula 
cerrada con otro compartimento vacío, y la válvula se abre para permitir el paso 
del gas, entonces éste se expandirá hasta ocupar todo el volumen disponible, y 
lo hará de forma espontánea. Además es un proceso irreversible: El proceso 
opuesto nunca sucederá de manera espontánea. 

A la hora de construir máquinas térmicas, que es lo que a un ingeniero le 
interesa de la termodinámica, han de tenerse en cuenta dos cosas: 

• De forma espontánea nunca puede existir un proceso que consista solo en 
flujo de calor desde un foco frío a uno caliente. 

• Tampoco sucederá un proceso que consista solo en conversión de calor en 
trabajo. 

5.7.2 Máquinas Térmicas 

Una máquina térmica es un dispositivo que funciona cíclicamente, con el propósito 
de convertir la máxima cantidad posible de calor en trabajo. En el desarrollo y 
estudio de máquinas térmicas se observó que no todos los procesos compatibles 
con el Primer Principio de la termodinámica son realizables, lo cual dio origen 
al Segundo Principio, que en definitiva es lo que limita el rendimiento de una 
máquina. 

El Segundo Principio de la Termodinámica, tal y como lo enunció Lord 
Kelvin: Es imposible que una máquina térmica funcione cíclicamente sin pro- 
ducir otro efecto que extraer calor de un solo foco realizando una cantidad de 
trabajo equivalente. 

A continuación se muestran dos esquemas de máquinas térmicas cuyo fun- 
cionamiento es imposible. 
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Como consecuencia del enunciado de Kelvin, sabemos que un ciclo sí que 
puede absorber calor de un foco (caliente) para realizar un trabajo. Pero no 
toda la energía absorbida se convertirá en trabajo. Las máquinas térmicas 
ceden parte del calor absorbido a un foco a baja temperatura, y el resto es lo 
que se convierte en trabajo. Esquemáticamente sería así: 
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5.7.3 Rendimiento en máquinas térmicas 

El rendimiento e de una máquina térmica a la que aportamos energía en forma 
de calor Qp, procedente de un foco caliente, para, cediendo calor Qp a un foco 
frío, obtener un trabajo M^ ne to aprovechable, se define: 

Wnetn 



Qc 



Es decir, la energía que obtenemos dividida por la que invertimos. 

En una máquina térmica que está funcionando cíclicamente, tras completar 
un ciclo se regresa a la situación inicial (de P, V , T, y también de energía 
interna U), con lo que Aí7 c j c i = 0. Por el primer principio: 



AU=Q=Q-W 



Q = W 



Observar que el balance de calor es el recibido (positivo) mas el cedido por 
el sistema (negativo): 

W = Q = Q c + Qf = \Qc\-\Qf\ 
Entonces, el rendimiento queda: 



W 



Q 



= i 



\Qi 



c 



\Qc 



5.7.4 El ciclo de Carnot 

El de Carnot es un ciclo reversible que realiza un gas ideal. Para regresar al 
estado inicial el gas atraviesa dos procesos isotermos y dos adiabáticos. El 
proceso se puede representar en un diagrama p-V como el de la figura. 




V 



Figure 5.6: Esquema del ciclo de Carnot 

El objetivo del ciclo de Carnot es la realización de trabajo hacia el exterior, 
mediante el aporte de energía a un ciclo. Parte de la energía que se aporta al 
ciclo en forma de calor absorbido procedente de un foco caliente, se transforma 
en trabajo que es aprovechable, mientras que el resto se convierte en calor que 
se cede a otro foco, éste de baja temperatura. 
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Teorema de Carnot: Ninguna máquina térmica que funcione cíclicamente 
entre dos focos térmicos puede temer un rendimiento mayor que la máquina de 
Carnot. La máquina de Carnot sigue un ciclo, llamado ciclo de Carnot, que 
consta de procesos de expansión y posterior compresión siguendo los siguientes 
pasos: 

1. Expansión isoterma y absorción de calor de un foco caliente. 

2. Expansión adiabática hasta alcanzar una temperatura inferior. 

3. Compresión isoterma y cesión de calor a un foco frío. 

4. Compresión adiabática para recuperar el estado inicial. 

Ahora que ya conocemos cómo es el trabajo realizado (o recibido) y el calor 
absorbido (o cedido) por un gas ideal en procesos isotermos y adiabáticos, pode- 
mos estudiar uno a uno todos los pasos que se llevan a cabo en un ciclo de Carnot. 
Denominamos Te a la temperatura del foco caliente y Tp a la del foco frío. 

Paso 1: Recordemos que en procesos isotermos, a temperatura Te, y pasando 
de un volumen inicial Va a uno final Vb'- 

Vf¡ 

Q 1 =W 1 = nRT c ln -f 

Va 

Por tanto, en primero proceso (entre los puntos A y B del diagrama p-V), 
el sistema absorbe calor (Q > 0) procedente del foco caliente y realiza trabajo 
(W = Q, W > 0) porque es una expansión, el volumen aumenta durante el 
proceso (Vb > Va)- 

Paso 2: El siguiente paso es la expansión adiabática, con un enfriamiento al 
pasar de Te a Tp < Te- En procesos adiabáticos, Q = 0, mientras que el 
trabajo es: 

W 2 = -nc v (T F - T c ) 

El trabajo es positivo, realizado por el sistema. Una expresión que después 
nos será útil, a la hora de calcular el rendimiento del ciclo, es la siguiente: En 
procesos adiabáticos, TV 7-1 = ele. Entonces, durante este proceso adiabático, 
en el que se pasa de Te a Tp en cuanto a temperaturas, y de Vb a Vp¡ en 
volúmenes, se tiene que: 

TeV^ 1 = TpV^- 1 

Paso 3: El ciclo continua con una compresión isoterma, proceso que tiene 
lugar a baja temperatura Tp, pasando de un volumen inicial Ve a uno final Vd. 
Entonces, el trabajo y el calor son negativos: 

Qz = W 3 = nRTp ln -£ 
Ve 

porque Ve > Vd y por tanto el logaritmo del cociente es menor que cero. 
El calor es negativo, el sistema está caliente y cede calor a un foco frío. El 
trabajo es negativo, el sistema recibe trabajo. 
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Paso 4: Para terminar el ciclo se realiza una compresión adiabática con calen- 
tamiento, pasando de una temperatura Tp a Te > Tp. Terminado este proceso 
el ciclo queda completado, regresa al punto A del diagrama p-V. Al ser un 
proceso adiabático, no hay transferencia de calor. El trabajo es: 

W 4 = -nc v (T c - T F ) 

Es una cantidad negativa, con lo que desde el entorno se está realizando un 
trabajo sobre el sistema. De manera análoga a la anterior adiabática, tenemos 
la relación entre volúmenes y temperaturas: 

TcVj- 1 = TpV^ 1 

Rendimiento del Ciclo de Carnot El rendimiento del ciclo es el trabajo 
total dividido por el calor que el sistema ha absorbido. El trabajo total es: 

W T = Wx + W 2 + W 3 + Wi 

Obsérvese que el trabajo en las dos fases adiabáticas es igual pero de signo 
opuesto, con lo que W 2 + W¿± = =>■ Wp = Wi + W 3 . 

Por otra parte el calor que el sistema absorbe es Q\, que ya calculamos en 
el paso primero. 

Pero además hemos visto que se puede calcular el rendimiento utilizando 
otra expresión: 

e = 1 _\Qil 

\Qc\ 

es decir, la unidad menos el cociente entre los valores absolutos del calor cedido 
al foco frío y el calor absorbido del foco caliente. En nuestro ciclo, el calor 
cedido el foco frío es Qp — Q 3 , mientras que el absorbido es Qc = Qi- 



*% 



nRTp ln |g 



\Qp\ = \Q 3 \=nRT F 

\Qc\=Qi = nRT c ln^ 

e=l ^ 

Ahora bien, como se vio en el estudio de procesos adiabáticos para gases 
ideales, TV 7 ^ 1 , con lo que se dedujeron dos expresiones para los procesos 2 y 
4: 



TcV^ 1 = TpV¿- x \ _ (VbV- 1 = (VcV- 1 Vb = Vc 

TcVJ- 1 ^TpV^- 1 } \V A ) \V D ) V A V D 

Es decir, el rendimiento del ciclo de Carnot es: 

e = 1 - -f- T en K 

Te 

Con lo cual, a mayor relación de temperaturas de los focos Tc/Tp, mayor 
será el rendimiento teórico de un ciclo de Carnot, pero siempre sin llegar al 
límite ideal de la unidad (el rendimiento nunca podrá ser uno, siempre es menor: 
e< 1). 



